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SECTION  1 


INTRODUCTION  AND  SUMMARY 

This  project  studies  the  atomic  structures  of  boundaries.  The  tool  used  in 
the  study  is  the  cluster  variation  method  (CVM)  of  statistical  mechanics,  with 
the  work  leading  to  information  on  excess  free  energy  (surface  tension),  excess 
entropy,  and  the  amount  of  adsorption  and  the  profile  of  the  density  across  the 
boundary.  This  project  is  unique  in  obtaining  such  information  from  a 
microscopic  point  of  view.  And  while  different  from  some  phenomenological 
approaches,  the  CVM  technique  is  thermodynamically  consistent,  in  the  sense  that 
thermodynamic  relations  are  always  satisfied. 

In  the  period  preceding  the  present  report  time  frame,  we  pursued  three 
different  directions.  The  first  considered  the  boundaries  in  binary  fee 
crystals,  the  second  examined  the  grain  boundary  structure,  and  the  third  was 
concerned  with  the  scalar-product  expression  of  the  boundary  free  energy.  In 
the  present  report  period,  the  work  in  the  first  two  categories  has  been  further 
exploited. 

(A)  In  working  with  the  statistical  mechanics  of  binary  fee  alloys  we  need 
a  tetrahedron  made  of  four  lattice  points.1*2  Our  two  previous 
papers3* 4  worked  on  the  antiphase  and  interphase  boundaries  (APB  and  IPB) 
in  the  CU3AU  ordered  phase  and  the  disordered  phase  using  a  tetrahedron.  We 
found  that  the  boundary  free  energy,  o  (of  both  APB  and  IPB),  goes  through  a 
maximum  when  plotted  against  the  temperature.  This  behavior  of  o  is  similar  to 
the  yield  strength  of  the  material  and  suggests  the  close  relation  between  the 
two  properties. 

The  work  on  the  fee  boundaries  has  been  continued  and  a  rotational  APB  in 
CuAu  ordered  phase  has  been  studied,  as  is  reported  in  Appendix  A.  Different 
from  our  previous  work  in  Refs.  3  and  4,  the  left-  and  right-hand  sides  of  this 
boundary  are  not  symmetrical.  Because  of  the  asymmetry,  the  boundary  behaves  in 
a  markedly  different  way  from  previous  boundaries.  By  contrast,  a  first  order 
phase  transition  occurs  within  the  boundary.  The  amount  of  adsorption  of  Au 
atoms  at  the  boundary  changes  sign  as  the  temperature  increases.  Another  kind 
of  phase  transition  is  also  derived  at  a  non-stoichiometric  composition. 


Although  the  theory  Is  thermodynamically  self-consistent,  we  decided, 
because  these  findings  are  so  unique,  to  spend  ample  time  to  verify  the  accuracy 
of  the  theory  before  publishing  the  result. 

(B)  In  previous  studies  of  the  structure  of  a  two-dimensional  grain 
boundary5  we  used  a  DSC  lattice  as  the  reference  lattice  on  which  atoms  can 
move.  The  lattice  constant  of  the  DSC  lattice  is  smaller  than  that  of  the 
crystals  on  both  sides  of  the  boundary.  We  found  that  the  grain  boundary 
changes  into  the  liquid-like  high  temperature  phase  at  the  temperature  range 
much  below  the  bulk  melting  point. 

We  have  extended  these  findings  and  have  worked  on  a  three-dimensional 
grain  boundary,  as  reported  in  Appendix  B.  The  work  is  an  ABC-CBA  stacking 
boundary  in  a  fee  crystal.  The  reference  lattice  has  three  times  as  many 
lattice  points  as  the  fee  lattice  of  concern.  The  basic  cluster  of  the  CVM 
method  is  a  three-layer  rhombus  prism  with  twelve  lattice  points.  Since  the 
structure  is  three-dimensional,  the  number  of  variables  is  large  (about  90,000). 
Nevertheless,  the  natural  iteration  method5  can  solve  the  set  of  simultaneous 
nonlinear  algebraic  equations.  For  computation  a  CRAY  computer,  which  is 
capable  of  vector  processing,  was  used. 

The  calculation  leads  to  thermodynamic  information  about  the  boundary:  the 
excess  free  energy,  the  entropy  and  the  density  profile.  Additional  work  is 
planned  for  different  temperatures  and  chemical  potentials,  and  for  (110)  and 
(112)  directions  of  the  boundary.  Thus  far,  however,  only  two  temperatures  for 
the  (110)  boundary  have  been  calculated,  and  further  commentary  on  the 
noteworthy  features  of  this  boundary  would  be  premature. 

(C)  In  Appendix  C  we  have  printed  a  still  unpublished  paper  by 

Dr.  John  W.  Cahn  of  the  National  Bureau  of  Standards  who  has  been  collaborating 
on  this  project.  This  paper  concerns  faceting  and  dissociation  of  grain 
boundaries  and  is  expected  to  be  an  important  work  in  the  field. 
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SECTION  2 


FIRST  ORDER  TRANSITIONS  WITHIN  ROTATIONAL  APB's  IN  Ll0  PHASE 

The  CuAu-type  ordered  phase  is  called  Ll0  and  is  a  layer  structure  with 
alternating  Cu  and  Au  layers*  We  placed  these  layers  perpendicular  to  the  plane 
of  a  paper  to  form  the  right-hand  side  of  the  boundary.  We  kept  the  layers 
perpendicular  to  the  paper  and  turned  them  90°,  then  placed  them  on  the  left- 
hand  side  of  the  boundary.  The  boundary  formed  is  illustrated  in  Figure  1  of 
Appendix  A. 

We  used  the  tetrahedron  as  the  basic  cluster  and  implemented  the  CVM2  to 
calculate  the  boundary  structure.  The  mathematical  technique  is  similar  to  the 
one  used  before4  in  calculating  the  APB  and  IPB  related  to  the  CujAu  phase.  The 
grand  potential  G  of  the  entire  system  including  the  boundary  region  is 
minimized,  keeping  the  chemical  potentials  fixed.  The  G  quantity  is  defined  as 

G  -  F  -  |PiNi  **  E  -  TS  -  (P1N1  +  P2N2)  ,  (1) 

where  F  is  the  Helmholtz  free  energy,  E  is  the  internal  energy,  S  is  the 
entropy,  y*  is  the  chemical  potential,  and  N^  is  the  total  number  of  the 
i-th  atom  (i»l  and  2).  In  minimizing  G,  we  placed  the  end  conditions  for  the 
far  left  and  far  right  ends  in  the  respective  homogeneous  state. 

In  the  formulation  we  did  not  include  vacancies  or  raany-body  interactions. 
Therefore,  the  results  are  symmetrical  with  respect  to  the  interchange  of  l  and 
2  (Cu  and  Au)  species.  When  we  express  ~  ^2  simply  as  y,  the  stoichiometric 
composition  (Cu  and  Au,  both  50%)  corresponds  to  U  ■  0.  We  may  choose  U  >  0  for 
the  Au-rich  side. 

When  we  calculated  the  excess  adsorption,  f,  of  Au  at  the  boundary,  we 
discovered  that  a  high  temperature  r  Is  positive  when  y  >  0  and  that  r  changes 
discontinuous ly  to  a  negative  value  when  y  crosses  y  ■  0  to  the  U  <  0  side. 

This  behavior  shows  that  there  are  two  phases  coexisting  at  y  *  0  for  this 
temperature  and  that  the  transition  is  first  order.  This  behavior  is  shown  in 
Figure  4(a)  of  Appendix  A. 
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For  the  low  temperature  region,  we  derived  the  first  order  transition  of  a 
similar  kind.  However,  the  low  temperature  behavior  is  distinctly  different 
from  the  high  temperature  one,  as  is  discussed  in  detail  in  Section  4  of 
Appendix  A;  we  can  define  a  transition  from  the  low  to  high  temperature  behavior 
at  kT3/|w|  ■  0.515  in  Figure  3  of  Appendix  A.  The  transition  at  T3  and  at  P  ■  0 
is  a  third  order. 

The  behavior  of  this  system  around  T3  satisfies  all  of  the  thermodynamic 
requirements.  However,  the  theoretical  results  are  so  unique  that  we  decided 
not  to  publish  the  findings  at  this  time,  but  rather  to  spend  more  time 
examining  the  validity  of  the  results. 

Basides  the  first  order  transition  discussed  above,  there  is  one  more  phase 
transition  near  the  composition  0.54,  as  shown  in  Figure  5  of  Appendix  A.  This 
phase  transition  is  due  to  the  order-disorder  change  within  one  of  the  Cu  layers 
of  the  layer  structure. 


SECTION  3 


THEORY  OF  ABC-CBA  STACKING  BOUNDARY  IN  fee  STRUCTURE 


The  (III)  planes  of  the  fee  structure  Is  stacked  as  ABCABC...  as  opposed  to 
the  ABAB...  stacking  of  the  hep  structure.  A  boundary  in  the  fee  structure  Is 
formed  when  the  stacking  is  ABC...  on  the  left-hand  side  and  CBA...  on  the 
right.  We  discuss  the  structure  of  this  stacking  boundary  in  Appendix  B. 

Atoms  were  not  only  allowed  to  be  located  on  fee  lattice  points,  but  also 
at  interstitial  positions.  The  interstitial  positions  are  defined  by 
translating  the  fee  lattice  perpendicular  to  the  (111)  direction,  as  shown  in 
Figure  2  of  Appendix  B.  Including  these  interstitial  positions,  the  reference 
lattice  forms  a  hexagonal  lattice  which  has  three  times  as  many  lattice  points 
as  the  original  fee  lattice.  This  hexagonal  lattice  plays  the  role  of  the 
reference  DSC  lattice  used  in  the  two-dimensional  grain  boundary  examinations 
discussed  in  Ref.  5. 

The  basic  cluster  is  the  three-layer  rhombus  prism  (defined  in  the 
reference  hexagonal  lattice),  as  shown  in  Figure  4  of  Appendix  B.  The 
subclusters,  V2,  U3,  U2,  P3  and  P2,  shown  in  Figure  4,  are  also  used. 

In  the  reference  hexagonal  lattice,  atomic  pairs  on  adjacent  lattice  pairs 
are  excluded.  The  second-neighbor  atomic  pair  is  the  stable  pair  and  is 
assigned  an  attractive  potential  energy,  -e  (e  being  positive).  The  potential 
energy  assignment  for  atomic  pairs  is  defined  in  Table  1  of  Appendix  B.  The  0-9 
pair  is  stable  in  the  fee  structure  and  the  0-8  pair  is  stable  in  the  hep 
structure.  In  view  of  our  work  on  the  fee  we  assigned  an  attraction  to  0-9  and 
a  repulsion  to  0-8. 

The  entropy  expression  based  on  the  rhombus  prism  is  one  of  the  key  points 
in  the  formulation.  It  is  derived  by  using  the  standard  derivation  of  the 
entropy  expression  in  the  CVM. 

The  equilibrium  state  is  derived  by  minimizing  the  grand  potential,  which 
is  of  the  form  of  Equation  (1),  the  difference  being  that  Nj  is  the  total  number 
of  atoms  in  a  system,  is  the  chemical  potential  of  the  atoms,  and  w2  *  0*  We 
impose  the  end  condition  where  the  left  end  of  the  system  is  in  the  homogeneous 
ABC  stacking  state  and  the  right  end  is  the  CBA  stacking  state. 
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Since  the  number  of  independent  variables  is  about  90,000,  it  is  necessary 
to  use  a  fast  computer.  A  CRAY  computer,  which  has  vectorization  capability, 
solved  the  variables  by  the  iteration  method  in  about  50  min  CPU  time. 

The  density  profile  across  the  boundary  is  shown  in  Figure  6  of  Appendix  B 
for  the  (110)  boundary  orientation. 

In  order  to  derive  useful  information  from  the  computation,  it  is  necessary 
to  work  out  several  different  temperatures  for  at  least  two  orientations  of  the 
boundary.  Further  work  along  this  line  is  being  planned. 


SECTION  4 


t: 


TRANSITIONS  AND  PHASE  EQUILIBRIA  AMONG  GRAIN  BOUNDARY  STRUCTURES 


N 


The  geometry  of  a  grain  boundary  is  characterized  by  three  angular 
variables,  R,  and  two  angular  variables,  N.  The  R  variables  specify  the 
misorientation  of  the  two  grains  meeting  at  the  boundary,  and  the  N  variables 
are  for  the  normal  to  the  boundary  plane.  The  excess  free  energy,  Y»  depends  on 
R  and  N.  Recently,  Dr.  John  W.  Cahn,  our  consultant  and  collaborator,  noted 
that  the  relation  between  Y  and  N  is  analogous  to  the  relation  between  the  bulk 
Helmholtz  free  energy  and  the  composition.  When  this  analogy  is  extended,  the 
well-established  concept  and  formulation  of  phase  equilibria  of  bulk  systems  can 
be  translated  into  the  new  concept  and  formulation  of  phase  and  phase  equilibria 
among  grain  boundaries.  In  this  way,  Cahn  has  developed  a  theory  of  grain 
boundary  faceting  and  dissociation.  This  is  a  pioneering  and  novel  idea  in  the 
field  of  grain  boundary  physics  and  will  undoubtedly  open  up  new  activities. 

Dr.  Cahn  kindly  agreed  to  have  this  unpublished  work  printed  in  this  report 
as  Appendix  C.  (The  work  was  done  at  the  National  Bureau  of  Standards. 

However,  the  research  in  Ref.  40  was  supported  by  this  contract,  as  acknowledged 
at  the  end  of  the  paper.) 
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APPENDIX  A 


First  order  transition  within  rotational  APB's 
in  Ll0  phase* 

Ryoichi  Kikuchi 

Hughes  Research  Laboratories 


ABSTRACT 


The  90°  rotational  antiphase  boundaries  (APB)  in  the  L1q  phase  (CuAu-type 
ordered  phase  in  the  fee  lattice)  are  studied  using  the  tetrahedron  approxima¬ 
tion  of  the  cluster  variation  method.  Calculation  is  done  for  a  nearest- 
neighbor  interaction  without  many-body  forces.  At  the  stoichiometric 
composition,  excess  adsorption  of  Au  at  the  boundary  Increases  discontlnuously 
as  the  boundary  changes  from  Cu-rich  to  Au-rich  lsothermally,  indicating  a 
first  order  transition.  There  is  another  first  order  transition  at  a  larger  Au 
composition.  The  nature  of  the  transitions  is  shown  by  plotting  profiles  of 
densities  on  four  sublattices. 


*  Supported  by  U.S.  Army  Research  Office 


A.  INTRODUCTION 

Properties  of  boundaries  In  the  fee  lattice  were  studied  before1  using  the 
tetrahedron  as  the  basic  cluster  In  the  cluster  variation  method  (CVM).2*3 
Calculations  were  done  on  the  antiphase  boundaries  (APB)  In  the  A3B  phase  (Ll2) 
and  on  the  interphase  boundaries  (IPB)  between  the  Ll2  and  the  disordered  (A3) 
phase.  The  excess  free  energies  due  to  the  IPB  and  APB  showed  a  marked  feature 
in  that  they  go  through  maxima  as  the  temperature  Increases. 

In  the  present  paper,  we  apply  the  same  tetrahedron  treatment  of  the  CVM  to 

APBs  in  the  AB-type  phase  (Ll0)  in  the  fee  lattice.  The  main  part  of  the  paper 

is  concerned  with  the  boundary  illustrated  in  Figure  1,  which  is  drawn  for  the 
Au-Cu  alloy.  The  Ll0  phase  has  a  layer  structure  of  Au-rich  and  Cu-rich  planes 
alternating  with  each  other.  A  line  in  Figure  1  indicates  a  crystal  plane  that 
is  perpendicular  to  the  plane  of  the  paper.  We  see  that  the  left  side  of  the 

boundary  is  rotated  by  90°  with  respect  to  the  right  side.  Since  the  layer 

planes  on  the  left  side  are  90s  to  the  direction  of  the  boundary,  we  may  call 
this  the  90°-rotational  boundary. 

As  shown  in  (A)  and  (C)  in  Figure  1,  two  cases  are  distinguishable, 
depending  on  whether  the  layer  closest  to  the  boundary  on  the  right-side  phase 
is  Au-rich,  as  in  (A),  or  Cu-rich,  as  in  (C). 

B.  GRAND  POTENTIAL  MINIMIZATION 

We  can  distinguish  four  sublattices.  We  number  the  lattice  planes  as 
shown  in  Figure  1.  The  odd  sublattices  I  and  III  are  on  odd  planes.  We  use  the 
index  1*1  for  a  Cu  atom  and  1*2  for  a  Au  atom.  The  probability  that  a  I  sub¬ 
lattice  point  at  position  n  is  occupied  by  the  i-th  species  is  written  as 
xi,n(i).  Since  vacancies  are  not  considered  in  the  formulation,  x's  are 
normalized  for  any  n  and  for  a  ■  I,  II,  III,  and  IV  as 
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Figure  1.  Structure  of  the  908-rotational  APB.  The  Au  and 
Cu  layers  are  perpendicular  to  the  plane  of  the 
paper.  In  (A),  the  Au  layer  on  the  right  is 
closest  to  the  boundary,  while  in  (C),  the  CU 
layer  is  closest  to  the  right. 


As  was  done  In  the  previous  work  on  the  fee  boundaries, 1  we  use  as  the 
basic  cluster  a  tetrahedron  made  of  four  adjacent  lattice  points.  There  is  one 
kind  of  tetrahedron  for  each  n.  Since  each  lattice  point  is  occupied  either  by 
Cu  or  Au,  there  are  16  possible  configurations  of  a  tetrahedron.  Probabilities 
of  these  configurations  are  used  as  the  basic  variables,  which  we  will  call  the 
z  variables. 

In  order  to  solve  the  equilibrium  structure  of  the  boundary,  we  first  write 
the  grand  potential  G  in  terms  of  the  z  variables,  and  then  minimize  G  with 
respect  to  z's. 

The  grand  potential  G  is  defined  as 

G  *  E  -  TS  ■  (y^  +  M2N2)  ,  (2) 

where  E  is  the  internal  energy,  S  is  the  entropy,  yj(i“l  and  2)  is  the 
chemical  potential  of  the  i-th  species  and  N  is  the  total  number  of  the  i-th 
atoms  in  the  system. 

The  energy  is  written  based  on  the  values  for  nearest  neighbor  pairs. 
Without  the  loss  of  generality  we  can  assign  the  energy  value  0  for  the  i-i  pair 
(i  ■  1  or  2)  and  w  for  the  1-2  pair,  as  was  done  in  Ref.l.  The  many-body 
interaction  parameters,  a  and  6,  used  in  Ref.l  however,  are  made  equal  to  zero 
in  the  present  paper  in  order  to  simplify  computations,  and  because  the 
Inclusion  of  a  and  B  does  not  change  the  qualitative  features  of  the  results. 

The  entropy  expression  is  the  key  element  of  the  CVM.  It  is  written  in 
terms  of  the  z  variables  and  of  the  probabilities  for  subclusters,  namely, 
lattice  pairs  and  lattice  points.  Triangles  are  subclusters  of  a  tetrahedron, 
but  do  not  appear  in  the  entropy  expression  of  the  tetrahedron  treatment.  There 
are  six  different  kinds  of  lattice  pairs  for  each  n,  depending  on  how  the  four 
sublattices  are  connected. 


In  writing  the  chemical  potential  terras  in  Equation  (2),  when  we  do  not 
include  vacancies  in  the  lattice,  the  sum  of  Uj  and  U2  vanishes  so  that  we 
define  the  independent  U  as 

=  y2  =  U  .  (3) 

Due  to  the  symmetry  of  the  problem,  M  *  0  at  the  50-50  composition,  and  U 
increases  with  the  Au  composition. 

The  equilibrium  structure  of  the  boundary  is  derived  by  minimizing  the 
grand  potential  G,  keeping  T  and  U  fixed,  with  respect  to  the  tetrahedron 
variables  under  the  normalization  constraints,  continuity  constraints  and  the 
end  structure  constraints.  The  continuity  constraints  guarantee  that  two  z's  on 
adjoining  position  share  a  lattice  pair.  For  the  left  end  conditions  we  require 
that  xj(i)  “  xjy(i)  and  xn(i)“  xjjj(i)  take  the  bulk  values  of 
the  Cu  and  Au  sublattices,  respectively.  For  the  right  end,  xj(i)“  xm(i) 
and  xn(i)  ■  xIV(i)  take  the  corresponding  values.  These  choices  are  in 
accordance  with  the  Figure  1  assignment. 

Differentiations  of  G  lead  to  simultaneous  non-linear  equations.  They  are 
solved  by  the  natural  iteration  method  (NIM),  which  has  been  developed  before.4 
Each  tetrahedron  has  16  z  variables.  Typically,  we  work  n  (the  position  of  the 
plane)  ranging  from  1  through  50,  leading  to  the  total  800  z  values.  The  con¬ 
tinuity  constraints  require  additional  150  Lagrange  multipliers.  In  spite  of 
the  relatively  large  number  of  variables,  the  NIM  converges  well.  The  Lagrange 
variables  are  determined  by  the  iteration  scheme  called  the  minor  Iteration. 1 

C.  BOUNDARY  PROFILE 

A  typical  profile  is  shown  in  Figure  2.  This  plots  the  occupation 
probability  x(2)  of  a  Au  atom  at  a  lattice  point.  The  boundary  is  around  the 
positions  n  *20  through  25. 


On  the  left  end  of  Figure  2,  Xjj(2)  and  Xjj^(2)  are  practically 
the  same  and  are  close  to  unity.  This  means  II  and  III  sublattices  are 
occupied  preferentially  by  Au  atoms  In  agreement  with  the  sublattice  assigment 
in  Figure  1.  As  we  come  closer  to  the  boundary,  xji(2)  and  x-£jj(2) 
gradually  deviate  from  each  other,  as  do  xj(2)  and  Xjy(2). 

The  structure  of  the  sublattice  occupation  on  the  right  side  of  the 
boundary  is  different.  At  and  beyond  the  position  n=22,  xj j ( 2 )  and 
xiy(2)  are  exactly  equal.  The  two  values  split  off  suddenly  at  n*20. 

Similarly,  xj(2)  and  xj^(2)  are  equal  at  and  beyond  n»23,  and  split  off 
at  n=21. 

We  now  see  the  difference  of  the  right  and  left  sides  of  the  boundary. 

Coming  from  the  right,  II  and  IV  split  suddenly,  while  on  the  left,  II  and  III 
gradually  deviate  from  each  other.  The  difference  is  understood  based  on  the 
fact  that  II  and  IV  are  on  the  same  plane  parallel  to  the  boundary,  while  II  and 
III  are  on  different  planes. 

Coming  from  the  right  in  Figure  2(a),  n=22  is  the  last  position  at  which 
the  Au  sublattices  II  and  IV  are  the  same.  Since  this  point  appears  often  in 
the  subsequent  discussions,  let  us  call  it  the  branching  point  and  write  it 
as  njjr(Au)  -  22.  From  the  splitting  of  I  and  III  sublattices  we  can  write 
that  njjr(Cu)  ■  23.  At  the  nj,r(Au)  position,  the  plane  is  a  Au-rich 
layer,  and  at  njjC(Cu)  the  plane  is  Cu-rich.  Since  n^r(Au)  <  njjr(Cu) 
in  Figure  2(a),  we  see  that  this  is  the  structure  of  Figure  1(A)  rather  than  1(C). 
Let  us  call  this  the  A-type  boundary,  implying  that  a  Au-layer  is  next  to  the 
boundary.  The  boundary  of  Figure  1(C)  will  be  called  the  C-type  boundary. 

T).  TRANSITION  FROM  A-  TO  C-TYPE  BOUNDARY 

The  boundary  structure  Is  either  the  A-type  or  C-type  and  not  in  between. 
Therefore,  we  expect  a  transition  between  the  two  types.  In  the  present  formu¬ 
lation  many-body  interaction  is  excluded  and  thus  the  system  is  symmetric  with 
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respect  to  the  Interchange  of  Cu  and  Au.  This  means  the  transition  from  A-  to 
C-type  boundary  occurs  at  the  50-50  composition  or  at  u  a  0. 

Figure  3  plots  the  excess  amount  of  Au  atoms  adsorbed  at  the  boundary 
region  in  the  A-type  boundary  at  the  50-50  composition.  The  value  of  r^u 
in  Figure  3  is  per  one  lattice  point  (for  the  area  of  a  )  in  a  plane  parallel 
to  the  boundary.  At  the  same  composition,  the  plot  of  excess  Au  atoms  in  the 
C-type  boundary  is  a  mirror  image  of  the  Figure  3  curve  with  respect  to  the 
center  vertical  line. 


Figure  4(a)  plots  schematically  the  excess  free  energy  a  of  the  boundary  at 
a  constant  temperature  against  P.  Note: 


3o 

3p 


r 

Au 


r 

Cu 


(4) 


Since  both  the  A-type  and  C-type  boundaries  exist  at  the  same  composition  and  at 
the  same  temperature,  and  since  3a/ 3p  changes  discontinuous ly  as  the  C-type 
boundary  changes  to  the  A-type  at  V  »  0,  the  transition  is  of  the  first  order. 

The  curve  in  Figure  3  crosses  the  center  vertical  line  at  about  kT3/|w|  ■ 
0.515.  At  this  temperature  the  remains  zero  when  we  go  from  the 
C-type  boundary  to  the  A-type  boundary,  and  hence  the  transition  is  acci¬ 
dentally  third  order. 

Below  T3  the  transition  is  again  first  order,  but  the  curves  behave 
differently,  as  shown  in  Figure  4(b). 

Figure  5  shows  another  demonstration  of  the  transition.  The  example 
is  at  kT/|w|  ■  0.7.  As  we  see  in  Figure  2(a)  and  (b),  the  Au  branching  point 
at  this  temperature  is  placed  at  n^r(Au)  ■  22.  Solid  curves  in  Figure  5 
are  x^(2)  and  xju(2)  at  n»21,  which  is  the  point  one  lattice  plane  to 
the  left  of  njjr(Au).  In  Figure  2(a),  x^(2)  and  xjh(2)  are  split 
at  n*21,  while  they  are  together  there  in  Figure  2(b).  Figure  5  shows  this 
behavior  by  the  splitting  of  xj(2)  and  xjjj(2)  at  x(2)-0.5. 
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Figure  3.  Adsorption  of  Au  atoms  at  the  boundary  plotted 
against  temperature.  The  curve  is  on  one 
branch  of  the  free  energy  surface  corresponding 
to  Figure  1  (A) . 
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Figure  4.  Schematic  diagram  of  o  versus  y  at 
a  constant  temperature.  (a)  is  for 
T  above  T^,  and  (b)  is  below. 
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Figure  5  indicates  another  phase  transition  at  about  x(2)*0.54.  For 
x(2)  <  0.54  the  branching  positions  are  nt,r(Au)  =  22  and  n^r(Cu)=23. 

For  x(2)  >  0.54,  the  latter  shifts  to  n^r(Cu)=25,  as  shown  in  Figure  6. 

It  is  worth  noting  that  I  and  III  sublattices  on  the  n=23  plane  are  split, 
although  the  plane  is  sandwiched  between  the  two  planes,  n=22  and  n*24,  on  which 
II  and  IV  sublattices  are  the  same. 

R.  DISCUSSION 

Figure  7  plots  the  excess  Au  atoms,  against  the  composition  at 

kT/|w|=0.7.  There  are  small  gaps  at  x(2)=i0.5,  0.542,  and  0.458,  and  they 
correspond  to  the  two  phase  transitions  discussed  in  Section  4  and  displayed  in 
Figure  5. 

Figure  8  plots  the  excess  free  energy  o  due  to  the  boundary  at  the  50-50 
composition  against  temperature.  This  curve  is  to  be  compared  with  Figures  6 
and  10  of  Ref.  1.  The  curves  in  Figure  8  are  raonotonically  increasing  as  T 
decreases,  since  this  is  for  a  constant  P(=0.0).  If  we  plot  a  for  a  constant 
composition,  other  than  x(2)*0.5,  we  expect  a  maximum  at  a  finite  temperature. 

For  general  information.  Figure  9  shows  \i  versus  x(2)  at  kT/|w[«0.7. 
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Figure  6.  Density  profile  of  Au  atoms  on  the  four  sublattices 
across  the  boundary  at  kT/|  w|  =  0.7  and  for 
U  =  1.3.  The  overall  density  of  Au  throughout  the 
system  is  0.549. 
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X(2),  DENSITY  OF  Au 


Figure  7.  The  excess  Au  adsorbed  at  the  surface  against  the  overall 
density  of  Au  throughout  the  system.  The  temperature 
is  kT/  w  =  0.70. 


Figure  8.  The  excess  free  energy  a  at  the 

stoichiometric  composition  plotted 
against  temperature.  Within  a 
crystal  plane  parallel  to  the 
boundary  the  area  per  lattice  point 
is  written  as  a2. 
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X(2),  DENSITY  OF  Au 


Figure  9 


The  chemical  potential  y  plotted  against 
the  density  of  Au  in  the  bulk  homogene¬ 
ous  phase. 
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THEORY  OF  ABC-CBA  STACKING  BOUNDARY  IN  fee  STRUCTURE* 

Ryolchi  Klkuchi 
Hughes  Research  Laboratories 

ABSTRACT 

In  the  fee  lattice  a  boundary  Is  formed  when  the  stacking  Is  ABC...  on  the 
left-hand  side  and  CBA...  on  the  right-hand  side.  The  structure  of  this 
boundary  Is  calculated  using  the  cluster  variation  method  (CVM).  A  hexagonal 
lattice  Is  used  as  the  reference  lattice;  the  original  fee  lattice  occupies  1/3 
of  the  reference  lattice.  A  three-layer  rhombus  prism  Is  chosen  as  the  basic 
cluster  of  the  CVM.  The  grand  potential  of  the  entire  system.  Including  the 
boundary  region.  Is  minimized  with  respect  to  the  basic  variables.  The  set  of 
simultaneous  nonlinear  algebraic  equations  for  90,000  variables  Is  solved 
Iteratively  using  a  CRAY  computer.  Properties  of  the  boundary.  Including  the 
excess  free  energy  and  the  density  profile  across  the  boundary,  are  obtained. 
The  calculation  done  so  far  is  for  one  temperature  and  for  the  (110) 
orientation.  Calculations  for  more  temperatures  and  for  (112)  orientation  are 
being  planned. 


*Supported  by  U.S.  Army  Research  Office 


1.  Introduction 


In  a  series  of  papers, 1-4  we  have  theoretically  studied  structures  of  phase 
boundaries*  The  theoretical  method  we  used  has  been  the  cluster  variation 
method  (CVM).5  One  of  the  most  recent  works  in  the  series  has  been  the  study  of 
a  two-dimensional  grain  boundary4  in  which  we  discovered  that  the  boundary 
changes  from  the  low-temperature  crystalline  state  into  a  high-temperature 
liquid-like  state  at  a  temperature  region  much  below  the  bulk  melting 
temperature.  In  the  study  we  used  a  square  DSC  lattice  as  the  reference 
lattice,  and  worked  with  a  model  where  atoms  in  the  stable  crystalline  structure 
occupy  one  out  of  five  DSC  lattice  points. 

This  paper  is  an  extension  of  Ref.  4  and  is  the  first  attempt  to  calculate 
the  three-dimensional  grain  boundary.  The  model  we  will  use  can  be  briefly 
summarized  as  follows.  The  fee  structure  has  the  ABCABC...  stacking  while  the 
hep  has  an  ABAB...  stacking.  In  the  fee  lattice,  when  the  left-hand  side  has 
the  ABC...  stacking  and  the  right-hand  side  the  CBACBA...  stacking,  a  boundary 
is  formed  between  the  two  sides.  The  structure  of  this  boundary  is  the  subject 
of  our  concern  here. 

Goodhew,  Tan,  and  Balluffi6  reported  an  experimental  study  of  this 
boundary.  We  will  supplement  these  experiments  and  will  elucidate  the  structure 
of  the  boundary. 


2.  Model 

The  (111)  plane  of  the  fee  lattice  is  a  two-dimensional  triangular 
sublattice,  which  is  called  B.  We  supplement  this  sublattice,  B,  by  two  more 
sublattices,  A  and  C,  which  are  linear  translations  of  B,  as  shown  in  Figure  1. 
We  call  the  (111)  plane  of  Figure  1  layer  1,  and  the  (111)  plane  above  it  layer 
2.  In  layer  2  the  lattice  points  are  directly  above  the  C  position  and  we 
supplement  them  by  A  and  B  triangular  sublattices,  as  shown  on  the  left-hand 
side  of  Figure  2.  The  A  lattice  points  in  layer  3  are  similarly  supplemented  by 
points  B  and  C.  Together,  lattice  points  A,  B,  and  C  form  the  reference  lattice 
on  which  atoms  can  be  located. 

In  Figure  2,  A,  B,  and  C  sublattices  are  marked  with  a  circle,  a  square, 
and  a  triangle,  respectively,  and  the  original  fee  lattice  points  are  indicated 
by  filled  marks.  The  left-hand  side  is  the  ABC...  stacking  and  the  right-hand 
side  is  the  CBA...  stacking.  The  center  of  the  boundary  is  indicated  by  broken 
lines. 
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Figure  1.  The  reference  lattice  points  A, 
B,  and  C  on  the  layer  "1"  of 
the  fee  lattice 


The  interaction  energies  are  assumed  to  be  as  indicated  in  Table  1.  When 
an  atom  is  at  0,  no  atoms  can  come  to  points  1  or  4.  The  distances  0-2  and  0-5 
are  the  same  and  are  the  nearest-neighbor  distance  of  the  original  fee  lattice. 

In  order  to  treat  the  ABC-CBA  boundary  problem  of  our  present  interest,  it 
is  necessary  to  include  the  third  layer  interaction  and  distinguish  the  ABC  and 
ABA  stackings.  For  this  purpose  we  include  the  0-8  and  0-9  interactions,  as  in 
Table  1.  The  0-9  interaction  is  favored  for  the  fee  stacking,  and  the  0-8  is 
favored  for  the  hep  stacking.  Since  our  system  is  fee,  we  choose  the  parameter 
>  0  (attraction)  and  ejj  <  0  (repulsion). 

In  formulating  the  mathematics  of  the  problem  we  number  the  lattice  planes 
parallel  to  the  boundary  as  n-1,  n,  n+1 , . . • ,  as  shown  in  Figure  2.  On  each 
plane  we  distinguish  nine  sublattice  points  *-l,...,9  as  shown  in  Figure  2.  The 
point  *“1  on  layer  1  is  marked  by  a  circle  and  <m2  by  a  square,  and  so  on.  Note 
that  vertical  to  the  layer  shown  in  Figure  2,  the  circle  points  form  a  line, 
while  the  square  points  form  a  line  of  their  own,  and  the  triangle  points  form  a 
line  of  their  own.  Thus,  for  example,  on  the  plane  n,  a  <«4  point  is  about  a 
*“1  point. 
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Table  1.  Definition  of  the  Interaction  Potentials 
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The  boundary  plane  in  Figure  2(a)  is  parallel  to  the  (110)  plane  of  the  fee 

lattice.  Another  boundary  shown  in  Figure  3  parallel  to  the  (112)  plane  can 

also  be  formulated  using  a  similar  technique  and  will  be  reported  later. 


3.  Variables 


In  treating  this  problem  using  the  CVM  we  have  chosen  a  three-layer  rhombus 
prism  as  the  basic  cluster.  It  is  true  that  a  three-layer  hexagon  prism  leads 
to  a  better  approximation,  but  it  was  not  used  because  it  would  need  larger 
memory  space  in  computers.  The  variables  we  use  in  the  formulation  are 
illustrated  in  Figure  4. 
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Figure  4.  The  shapes  of  the  clusters  used  In  the  CVM  treatment. 

a.  V2  and  V3  Variables.  Configurations  of  a  rhombus  are  numbered  as  In 
Table  2.  Note  that  we  distinguish  three  orientations  of  the  rhombus:  left- 
tilted,  horizontal,  and  rlght-tllted.  Using  the  numbers  in  Table  2,  we  write 
the  probability  variable  of  a  three-layer  rhombus  prism  as  V3n> 

The  subscripts  1,  j,  and  k  take  one  of  the  numbers  in  Table  2,  1  being  for  the 
bottom  layer  configuration,  and  j  and  k  for  the  middle  and  top  layers, 
respectively.  A  rhombus  has  two  diagonals.  The  two  points  at  the  ends  of  the 
longer  diagonal  lie  on  the  same  k  lattice  points  in  Figure  2.  The  subscript  k 
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Table  2.  The  Definitions  of  the  Rhombus  Configuration 
Numbers.  A  dot  represents  an  atom. 
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1 


n 


is  this  <  number  for  the  end  points  of  the  longer  diagonal.  Because  the 
orientations  of  the  three  rhombi  (on  the  three  layers)  of  a  prism  are  the  same, 
we  have  a  restriction  on  the  i,j,  and  k  that  all  three  numbers  are  in  one  of  the 
three  columns  in  Table  2;  numbers  of  different  columns  of  Table  2  do  not  mix. 

Because  of  the  exclusion  requirement  in  Table  1,  when  the  bottom  layer  has 
configuration  2  (in  Table  2)  the  middle  layer  cannot  be  2  or  6.  This  exclusion 
relation  is  written  in  Table  3. 


Table  3.  The  Exclusion  Matrix 
MXC1.V  (I,J) 
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In  addition  to  the  V3  variables  we  use  V2n>jc.j>j  for  the  2-layer 
rhombus  prism.  The  meaning  of  the  subscripts  are  similar  to  those  for  V3. 

The  precise  definition  of  the  V3  variable  is  that  out  of  N  rhombi  of  the 
(n,k)  kind,  the  number  of  rhombi  which  have  the  configuration  (i,j,k)  is 
v3n,k;i, j,k* 

The  V2  variables  can  be  evaluated  when  the  V3  variables  are  given  as 


V2 


n,k;i,  j 


X 

k 


V3 


n,k;i, j,k 


.£  V3 

k  n,kM3;k,i,j 


(3.1) 
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where  kM3  Is  the  k-number  of  the  lattice  point  one  layer  below  k,  and  is 
calculated  as 

kM3  -  MOD(k  +  5,9)  +  1  .  (3.2) 

b.  V2  and  V3  Variables.  There  are  two  triangle  orientations,  left¬ 
pointing  and  right-pointing.  The  configurations  of  a  triangle  are  numbered  in 
Table  4.  The  probability  variables  for  the  three-layer  triangle  prism  are 
written  as  i , j ,k»  where  k  is  for  the  left  or  tight  tip  of  the 

triangle.  The  subscripts  i,j,  and  k  are  for  the  configuration  of  Table  4,  and 
no  numbers  on  different  columns  of  Table  4  mix.  The  exclusion  requirements 
corresponding  to  Table  3  are 

i  j*  j  and  j  j4  k  except  i  *  j  ■  1  and  j  *  k  s  1  .  (3.3) 

Table  4.  The  Definition  of  the  Triangle 
Configuration  Numbers.  A  dot 
represents  an  atom. 
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For  the  two-layer  triangle  prism,  we  define  V2n  ^.^j.  They  are 
derived  from  V3's  by 


V2 


n,k;l,  j 


Xv3 

k  V  n,k;t,j,k 


v3 

k  n,kM3;k,l,j 


(3.4) 


where  kH3  Is  defined  In  (3.2). 


c.  P2  and  P3  Variables.  These  variables  are  defined  In  Table  5.  The 
exclusion  requirements  are  the  same  as  In  (3.3).  The  reduction  relations 
similar  to  (3.3)  hold. 


Table  5.  The  Definition  of  the  Pair  and  Triplet 
Variables.  A  dot  is  an  atom. 
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d.  Reduction  of  U3  from  V3.  The  variable  V3n>k.  j>k  can  be  derived 
from  either  of  the  three  orientations  of  V3.  As  a  preparation  of  this  reduction 
we  write  the  relations  among  two-dimensional  triangles  (written  as 
in  Table  4  and  two-dimensional  rhorabl  (written  as  yn,k;i)  *n  Table  2. 

They  are 

xn,k;*=  [XLYLU,i  yn-l,kp2;i 

"  g  1M]t;iyn,k;W  . 4  (3*5a) 

g  [XLYR]A;1  yn>kpl;1+12 

for  the  left-pointing  triangle,  and 

Xn,ki«M  *  ^  (“"•It-i  y„,kpi;i 

*  S  [XRICI*u  Vnw  f°r  *-1 . 4  <5-5b> 

6 

£  [XRYR] ±  yn+1>kp2;i+i2 

for  the  right-pointing  triangle.  The  matrices  [  are  defined  in 

Table  6.  The  subscript  kPl  and  kP2  are  defined  as 
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Table  6.  The  Definition  of  the  Matrices  to  Connect  Rhombus 
Configurations  (l)  to  Triangle  Configurations  (&) 
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1 


kPl  *  k  +  1 


except  kPl  *  k  -  2  when  M0D(k,3)  ■  0  , 


(3.6a) 


and 


kP2  =*  k  -  l  , 


except  kP2  =  k  +  2  when  M0D(k,3)  *  1 


(3.6b) 


Using  the  matrices  defined  in  Table  6,  we  can  write  the  reduction  relations 
of  U3's  from  V3's  as 


o 

U3  ,  .  -  Y*  [XI.YL3]  .  .  ,  .  V3 

n,k;i,m,n  £  £w=1  J*,m,n;i,j,k  n-l,kP2;i, j,k 


i»S-i  txLYC3lz»m»n;i*j*kV3n»k;i+6»j+6»k+6 


(3.7) 


i.gui  lXLYR3]^,m,n;i,j,k  V3n,kPl;i+12,j+12,k+12  » 


where  we  define 


[XLlnJlt,m,.il,J,k  *  •  <3'8) 


and  similar  relations  for  [XLYC3]  and  [XLYR3].  In  the  subscripts  of  U3  on  the 
left-hand  side  of  Equation  (3.7),  n  appears  twice,  first  for  the  lattice 
location  and  then  for  the  triangle  index  (n*l,...,4),  but  no  confusion  is 
expected.  In  Equation  (3.7)  £,m  and  n  are  from  1  through  4. 
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The  relations  corresponding  to  (3.5b)  are 


U3n,k;*+4,m+4,n+4  *  t  *,m,n;i,  j,k  V3n,kPl;l,  j,k 


(3.9) 


and  two  more  expressions. 


e*  Reduction  of  P3  from  V3.  The  variables,  P3's,  can  be  written  as  linear 
combinations  of  V3's.  For  example. 


P3 


n,k;  £,m,n 


6 

V  [PV2] .  .  ,  ,  V3  .  .  .  , 

l,jTk-l  n,k;i,j.k 


(3.10) 


This  Is  for  the  top  lattice  point  where  the  atom  sits  In  configuration  2  on  the 
left  column  of  Table  2.  We  can  write  similar  relations  for  configurations  2,  3, 
and  4  of  Table  2  using  matrices  [PV3],  [PV4],  and  [PV5].  These  matrices  are 
defined  In  Table  7. 


4.  Grand  Potential 

The  mathematical  procedure  of  calculating  the  equilibrium  state  is  to 
minimize  the  grand  potential,  G,  keeping  the  chemical  potential  fixed.  The 
function  G  is  defined  as 


G  -  B  -  TS  -  HNa  ,  (4.1) 

where  M  Is  the  chemical  potential  and  is  the  total  number  of  atoms  in  a 
system. 


a*  Energy  B.  Calculation  of  the  energy  is  defined  in  detail  in  the  model 
in  Section  2.  Using  the  V3  variables,  the  energy,  E,  is  written  as  follows: 


9  6 

E-NEE  E  Sik  E  V3n  k-i+m  14m  k-Hn 

nk-1  1,J,k  m-0,6,12  n,k,i+m, j+m,k+m  , 


(4.2) 


where  N  Is  the  number  of  three-layer  units  in  a  system. 
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The  parameter,  Gt » j » It »  *s  t^ie  energy  per  rhombus  and  is  made  of  three 
parts: 


«<*>  -  +  e<l>  +  e<2>  +  e<2>  +  c<3> 

i  J  k  ij  jk  ik 


(4.3) 


where  Is  the  energy  within  a  layer,  ej2^  is  the  energy  between  adjacent 

(3) 

layers,  and  r-^  is  the  energy  between  the  first  and  third  layers. 

The  energy  due  to  the  diagonal  pair  of  a  rhombus  (for  the  i*6  case  in  Table 
2)  is  counted  in  three  V3's  (appearing  in  three  different  layers  of  a  V3). 

Thus 


ev 

i 


1) 


(4.4) 


A  pair-wise  Interaction,  e,  between  adjacent  layers  (for  the  pair  0-5  in 
Table  1)  is  counted  ten  times.  Thus 


1_ 

10 


(ED) 


U 


where  the  matrix,  (ED)^j,  is  defined  in  Table  8. 


(4.5) 


(3) 

The  third  layer  Interaction,  ,  is  due  to  either  “GEF  or  -eejj  in 

Table  1.  The  vertical  bond,  ~EE„>  is  counted  12  times:  4  times  in  the  rhombus 

H 

of  one  of  the  three  orientations  (Table  2).  The  fcc-type  interaction,  -eeF, 


is  counted  5  times,  i.e., 


a  half  of  e 


(2) 

ij  * 


Thus,  we  write 


c<3> 

ik 


12  MH  <EH>ik  +  5  «F  (ED)ik 


(4.6) 


where  both  matrices  are  defined  in  Table  8 


b.  Entropy  S.  The  CVM  entropy  expression  based  on  the  three-layer  rhombus 
prism  can  be  derived  from  the  entropy  expression  for  the  two-dimensional 
triangular  lattice  based  on  a  rhombus  cluster.  The  latter  was  derived 
previously  1  as 

S/k  (Pair)2 _ 

e  {Rhombus)  (Point) 


{Triangle) 
{Rhombus)  {Pair 


(4.7) 


where  c  is  an  integer.  When  the  diagonal  pair  of  the  rhombus  does  not 
contribute  to  the  potential  energy,  the  ratio  Inside  [  ]c  becomes  unity. 

When  the  diagonal  pair  does  contribute,  as  in  the  present  case,  we  should  use 
c  *  2.  Thus ,  in  the  present  case,  the  basic  relation  is 


_ (Triangle ) _ 

2 

{Rhombus)  {Point) 


(4.8) 


The  entropy  expression  for  the  three-layer  rhombus  prism  is  then  derived  as 


{3-ft.tr.pr.  )^ _ t _ {2-ft.tr.pr.  )^ _ 

3-^.rh.pr.  j  {3- £. triplet )  {2-&,rh.pr.)  (2-£. doublet . 


where  we  abbreviated  *,  «  layer,  rh.  =  rhombus,  pr 


prism,  and  tr.  =  triangle. 


Writing  in  full,  we  obtain 


S  =  kN 


9  f  6  2 

EE  -  E  X>(v3nk- 

n  k=l  i,j,k=l  m=0  *  ’ 


i+6ra, j+6m,k+6ra 


) 


+  S-l  S^^n.kji+bm.jW 
+  2  i,S.  £  ^(U3n»k;i+4m,j+4m,k+4m) 


4  1 


2  i5l  5  ^(U2n,k;i+4m,j+4m) 

2  2 

-  £  k-i  1  k>  +  £  «P2„  k-i 

i,j,k=l  n,k,l,  j,k  1>jH>  n.k.J 


where  we  define 


#(X)  =  X&nX  -  X 


(4.10 


(4.11 


In  Equation  (4.10),  N  is  the  number  of  three-layer  units  in  a  system  as  it  was 
in  Equation  (4.2). 


% 


N 


c.  Chemical  Potential  Term.  In  the  chemical  potential  term,  -PN^  in 
(4.1),  we  need  the  number  of  atoms,  N^,  which  we  write  as  follows.  Each  atom 
Is  counted  12  times  in  a  3-layer  rhombus  prism  of  one  orientation  (left-tilted, 
horizontal,  or  right-tilted,  as  in  Table  1).  Therefore, 


N  -  N 
A 


n,k; i+6m, j+6m,k+6m 


(4.12) 


where  Aj  is  the  number  of  atoms  in  the  rhombus  configuration  (i  *  1,...,  6, 
in  Table  2),  and  is  defined  as 


*>  0  for  i  =  1 

Aj  ■  1  for  i  =  2,  3,  4  and  5 


■  2  for  i  =  6 


(4.13) 


5.  Basic  Equations 

The  equilibrium  state  of  the  boundary  is  derived  by  minimizing  the  grand 
potential,  G,  and  keeping  the  chemical  potential,  W,  fixed.  In  so  doing, 
several  comments  are  in  order. 

a.  Choice  of  Independent  Variables.  We  choose  the  V3  variables  for  the 
basic  cluster  as  Independent.  Other  variables,  V2,  113,  U2,  P3  and  P2,  are  all 
for  the  subclusters  of  the  basic  cluster  and  can  be  written  as  linear 
combinations  of  the  V3  variables,  as  was  shown  in  Section  3. 

b.  Constraints  on  the  Independent  Variables.  Although  we  choose  V3's  as 
the  independent  variables  In  minimizing  G,  not  all  V3's  are  actually 
independent.  The  mutual  relations  to  be  satisfied  among  V3's  are  of  three 
kinds.  The  example  of  the  first  kind  is  the  equality  in  Equation  (3.1).  The 
constraints  of  the  second  kind  come  from  the  three  ways  of  writing  V3,  as  in 
Equation  (3.7),  and  additional  three  ways  In  Equation  (3.9).  The  third-kind  of 
constraint  is  the  normalization  relations: 
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-  n,k;  i+6mf  j+6m,k+6m 
» J  >  k3* 

for  each  n,k,  and  m30,  1  and  2 


(5.1) 


In  minimizing  G  with  respect  to  V3's,  we  use  the  Lagrange  multipliers,  a' s 
Y's  and  X's,  for  the  three  kinds  of  constraints. 


c. 

follows 


The  Basic  Equations.  The  basic  equations  for  V3s  are  written  as 
(we  use  i, j,k=l , . . . ,6 . ) : 


V3n,k;i,j,k  ~  V3n,k;i,j,k  eXp 


(5.2) 


where  the  exp  (  )  factor  is  determined  from  the  normalization  relations  in 

A 

Equation  (5.1).  The  un-norraallzed  part,  V3,  is  written  in  two  steps.  In  the 
first  step  we  include  the  Lagrange  multipliers,  a's  and  Y's,  mentioned  in  b. 
above,  and  write 


inV3  ,  ,  ,  .  -  £nV3  ,  ...  -  a  .  .  .  +  a 

n,k; i, j,k  n,k;i,j,k  n,k;i,j  n,kP3;j,k 


i,ra,n=l 


|  YLLn-4-l,km2;  i,m,n  [*L^k3]  i,m,n;  i,  j,k 


+  YRLn,kml; i+4 , m+4 , n+4 1 XLYR3 J  £,m,n;i, j,k  ( *  (5, 


The  V3  part  does  not  include  Lagrange  multipliers,  and  is  written  as 


*nV3n,k;l, j,k  “  3ei,j,k  +  ep(Ai  +  Aj  +  Ak)/36 


(5.4) 


+  -7  inV  ,  ,  ,  ,  +  inU  .  ...  -  7  inP  .  .  , 

2  n,k;i,j,k  n,k;i,j,k  4  n,k;i,j,k 


(5.4) 


The  j  ^  term  comes  from  the  energy  expression,  and  the  A  terras  come 

from  the  chemical  potential  term;  we  define  0  »  1/kT.  The  last  three  in  terras 


n,k;i, j,k 


~U  I  V  £. 


y  n,k;i,j  n,kP3;j,kJ 


AnUn  k-i  1  k  =  £ 

n»k, i, j,k  *m>n=1 


[XLYL3K  .  ,  .*n 
*,m,n;i,  j,k 


U3 


n.km2: &.m.n 


^2n+l,km2;  £,ra^2n+l,km2P3;m,n 


U3 


n .  kml ;  &+4 „  nr+4 .  n+4 


^2n,kml,  £+4,m+4^2n,kmlP3;nrt-4,n+4  ^ 


1/2 


E  ,  CPV2] 
in,  n=l 


P3 


in 


o.ki&.m.n 


^,in,n;i,j,k  .1/2 

1  n,k;  A,ra  n,kP3;m,nJ 


+  Ipv3]j  ,  ,  , 

JA,ra,n;l, j,k 


P3 


In 


n+l.k: fl.m.n 


^P2n+l,k;  £,m  P2n+l,kP3;m,n^ 


1/2 


+  [PV4^,m,n;i,j,kAn 


P3 


n.kPl;  E.m.n 


[P2n,kPl;*,m  ^n.kPlPSjm.n1 


1/2 


+  [PV5] »  .in 

J*,ra,n;i,  j,k 


P3 


a+l.JKP.2;  ^.a.n 


[P2n+l,kP2; i,m  P2n+l,kP2P3;m,n1 


1/2 


C 


We  derive  similar  expressions  for  V3n f ^ f j+£ f ^5  and 
v3n,k;i+12, j+12,k+12»  but  tlie  exact  equations  will  not  be  listed  here.  In 
the  latter  we  use  YLR  and  YRR  corresponding  to  Y's  in  Equation  (5.3).  The 
expression  for  v3n,k;i+6, j+6,m+6  contains  all  four  kinds  of  Y. 


6. 


Iterative  Solution 


We  solve  V3's  from  the  basic  equations  and  constraints  using  an  iterative 

o 

procedure  called  the  natural  Iteration  method  (NIM).  The  NIM  is  separated  into 
the  major  iteration  and  the  minor  iteration. 

a.  The  Major  Iteration.  (1)  We  start  with  the  input  values,  V3’s. 

(2)  The  subcluster  variables  V2,  U3,  U2,  P3  and  P2  are  derived  from  V3  by  the 
reduction  procedure  of  Section  3.  (3)  Using  (5.4)  and  (5.5),  we  calculate  V3. 

A 

(4)  The  un-no cmallzed  V3  in  (5.3)  is  written  using  the  Lagrange  multipliers  a's 
and  Y's.  Using  the  minor  Iteration  process  to  be  explained  in  (b)  below,  the 

A 

values  of  a's  and  Y's  are  determined  so  that  the  constraints  on  V3's  are 

satisfied.  (5)  After  the  minor  iteration,  the  normalization  factor,  X's,  are 

determined  from  Equations  (5.2)  and  (5.1)  to  end  up  with  the  output,  V3,  of  one 

major  Iteration  step.  This  output  is  used  as  the  input,  V3,  in  step  (1)  above. 

b.  The  Minor  Iteration.  We  determine  a's  and  Y's  in  Equation  (5.3)  for 
the  given  V3  in  such  a  way  that  the  first  and  second  kinds  of  constraints 
mentioned  in  Section  5b.  are  satisfied.  When  we  use  Equation  (5.2)  in  the 
first-kind  constraint  (3.1),  we  write 


exp 


^  an,k;  i+6m,  j+6ra^  exP^°n,k*t  i+6m,  j+6m^^n,k;  i+6m,  j+6m,k+6m 


exp 


^n.k;  i+6ra,  j+6m^  exP^  n,k;  i+6m,  j+6m^^n,km3; k+6m,  i+6m,  j+6m 


m  =  0,  1  or  2 


(6.1) 


Because  V3's  satisfy  the  normalization  of  Equation  (5.1),  we  can  choose  a  for 
i=j^i  a8  zero  without  the  loss  of  generality: 


a 


n,k; l+6m, l+6m 


0 


m  ■  0,  1  or  2 


(6.2) 
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When  we  use  Equation  (6.2),  we  can  eliminate  the  normalization  factor,  A,  by 
forming  a  ratio  of  Equation  (6.1).  We  define 


_iL  - 

«.  =  V  V3 

n,k; i+6m, j+6m  n,k; i+6m, j+6m,k+6m 

A  - 

*n,k; i+6m, j+6m  "  V^n,kM3;k+6ra,i+6m, j+6m  *  (6.3) 


Then  we  solve  Equation  (6.1)  in  the  iterative  form  as 


a(out)  ,  a(in) 

n,k; i+6m, j+6m  n,k; i+6m, j+6m 


l_k  rP^n.ki  i+6m.  1+6m  P^n.k:l+6m.l+6m~l 
2  n  [_®^n,k;  l+6m,  l+6m  PPn,k;  i+6m,  j+6m  J 


(6.4) 


where  ba  is  a  damping  factor.  When  ba  ■  0,  a(°ut)  is  equal  to 
a(ln).  ba  -  1  is  the  case  of  no  damping.  In  the  actual  computation  we 
used  b0  *  1/4. 

Using  a  procedure  similar  to  the  above  for  a,  we  derive  the  minor  iteration 
relations  for  Y's  as  follows!  In  order  to  satisfy  the  two  constraint  equations 
in  (3.7),  we  derive 


Yu/0^}  -  YLL 

n,k; *,ra,n 


f .  \  -  ("glUj  -  .  GLLR  -  *  .  > "] 

(in)  +  1  h  t  __  n,kj  i,m,_n  _  n.k;  1.1.1 

2  T  |au.n,kil.l.l  GL“n,ki«,m.nJ  •  <6‘: 


and 


YLR(°“CJ  -  YLR 
n,k;  *,ra,n 


where  we  define 


* .  \  i  GLRL  .  a  GLRR  1111 

(in)  ,  L  .  a,. _ n,vjA.JUja _ nJs*l*l*JL 

n,k;*,m,n  2  Y  GLRL  .  ,  ,  ,  GLRR  ,  .  » 

*  ’  *  n,k;l,l,l  n,k;*,ra,n 


(6.5) 


GLLL  a  n  =  V  [XLYL3] .  ,  ,  .  V3  .  .  .  .  , 

n,k;*,m,n  .  —j'  .  1  *,ra,n;i,j,k  n-1 ,kP2; i, j,k 

i »  3 » 

6 

GLLRn,k;*-,m,n  "  l  lXLYC3l  4,ra,n;  i,  j,k  V3n,k; i+6,  j+6,k+6 

6 

GLRLn,k;*,m,n  =  ±  [XLYR3]  A,m,n;  i,  j,k  V3n,kPl;  i+12,  j+12,k+!2 

GLRRn,k;*,m,n  =  t  5^  IXLYC31  4>m,n;  i,  j,k  Y3n,k; 1+6, j+6,k+6  *  (6,6) 

The  coefficient,  by,  in  Equation  (6.5)  is  a  damping  factor  of  similar  nature 
as  ba  la  Equation  (6.4).  In  the  actual  calculations  we  used  by  *  1/2. 

We  have  two  more  constraint  relations  for  V3  in  the  right-pointing  triangle 
variables  in  Equation  (3.9).  These  two  constraints  derive  YRL  and  YRR,  but  the 


details  are  omitted. 
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7 .  Homogeneous  Phase 


Before  calculating  the  boundary  structure,  we  solve  the  equilibrium  state 
of  the  homogeneous  phase.  We  can  use  the  same  equations  designed  for  the 
boundary  and  simply  make  the  variables  independent  of  the  lattice  location,  n. 

The  homogeneous  system  is  made  up  of  three  sublattices.  Referring  to 
Figure  2,  the  three  are 


\  "  1.  5,  9, 

<2  =  2,  6,  7, 

<3  -  3,  4,  8,  .  (7.1) 

When  the  temperature  is  very  high,  the  three  sublattices  become  equivalent  and 
the  system  is  disordered.  For  low  temperatures  when  the  sublattices  are  not 
equivalent,  we  can  define  the  long-range  order  as  the  difference  of  the  atomic 
densities  on  the  lattices  <2  and 

5  s  P(<2)  -  PC^)  .  (7.2) 

Note  that  we  make  <2  the  sublattice  on  which  atoms  lie  preferentially  as  in 
Figure  2.  Figure  5  shows  how  the  long-range  order,  5,  changes  as  a  function  of 
temperature.  We  chose  the  chemical  potential  as  u/e  «  -40. 

A  typical  calculation  of  kT/e  ■  4  and  u/e  =  -4  was  timed  as  follows.  The 
GPU  time  used  by  the  VAX/VMC  computer  was  200  minutes  for  450  major  iterations. 
The  corresponding  calculation  was  done  in  ten  minutes  using  the  CRAY  computer  of 
U.S.  Air  Force  Weapons  Laboratory. 
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Figure  5.  The  long-range  order  parameter  in  the  homogeneous 
bulk  phase  plotted  against  the  temperature  at 
M/e  -  -4.0. 

8.  Boundary  Calculations 

We  use  the  results  of  the  homogeneous  phase  calculations  as  the  end 
condition  of  the  boundary  calculations*  We  choose 

kT/e  -  4  and  u/e  -  -4  ,  (8.1) 

which  are  for  the  long-range  order  of  0.85,  as  we  see  in  Figure  5. 
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a.  Width  of  the  System.  For  the  boundary  computations  we  choose  the 
lattice  plane  number,  n.  In  Figure  2  to  go  from  l  through  ny,  the  center  being 
at  nc  =  (nw  +  l)/2.  If  nw  is  too  small,  the  end  effect  of  the  system 

makes  the  calculated  results  unreliable.  If  nw  is  too  large,  the  memory  space 
of  the  computer  cannot  contain  a  LI  of  the  variables  we  need.  As  a  compromise  we 
chose 

n  =17  and  n  =  9  .  (3.2) 

w  c 

b.  Left-Right  Symmetry.  In  the  actual  computations,  we  do  not  need  to 
work  for  n  from  1  through  riy.  Since  the  left-  and  right-hand  sides  of  the 
system  are  symmetric,  we  can  work  for  n  from  1  through  nc  +  l  and  use  the 
left-right  symmetry  relations.  We  note  the  following  kinds  of  symmetry: 

vl  =  VI 

9-n, <v;k+( 1 *1)6, j+( 1 ±1)6, i+( 1  ±1)6  9+n, <; l+(l +1)6, j+( 1 +l)6,k+(l  + 1)6 

for  i,j,k  =  l,...,  6  .  (8.3a) 

These  are  for  the  left  and  right  columns  of  Table  2.  For  the  center  column  of 
Table  2,  we  have 


V3  =*  V3  « 

9-n, <v;k'+6, j’+6,i'+6  9+n, i+6, j+6,k+6 

In  these  expressions,  <V  is  a  function  of  <  and  takes  the  value 


(8.3b) 


«V  =  4,5,6, l, 2, 3, 7, 8, 9  ,  (8.4) 

corresponding  to  <  8  1,....,9.  The  primes  on  i',  etc.  in  Equation  (8.3b)  are 
defined  as 


except  i'  ■  2  when  i  *  3 


and 


i'  =  3  when  i  ■  2 


(8.5) 
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The  symmetry  relations  for  a's  are 


N 


S 


V. 

»_  • 

ni 


tm 


9-n,  «*; j+(2-m)6, i+(2-m)6  a9+n,k; i+6m, j+6m 

for  i,j  *  and  m  «  0,1,2  .  (8.6a) 

The  subscript,  <<*,  is  defined  as 

ka  -  7, 8, 9, 4, 5, 6, 1,2,3  ,  (8.6b) 

corresponding  to  <  *  l, '...,9.  Note  the  minus  sign  on  a  on  the  right-hand  side 
of  (8.6a). 

The  symmetry  relations  for  Y's  are 

YLL_  ,  *  YRR_  n 

9-n,<v;n,m, *  9+n,  *,m,n 

YRL  -  YLR 

9-n,Kv;n,m, *  9+n,*;l,m,n  (8.7) 

where  we  use  <v,  defined  in  Equation  (8. A). 

c.  Initial  Values.  To  start  the  iteration,  the  initial  values  of  the  V3 

variables  are  chosen  as  follows.  For  n  -  1,...,  nc  -  1,  V3’s  are  set  equal  to 
the  homogeneous  phase  values.  At  the  center,  n^.,  we  choose  V3  as  the  average 
of  the  left-side  homogeneous  phase  value  and  its  symmetry  value  (calculated  by 

the  relations  in  Equation  (8.3).  We  use  the  same  idea  in  choosing  the  initial 

values  of  a's  and  Y's. 

d.  End  Condition.  After  the  iteration  has  started,  we  choose  the  end 
condition  at  n  ■  1  and  2  so  that  the  values  of  V3  at  n  -  1  and  2  are  the  same. 

e.  Computer  Time  and  Memory.  When  we  choose  r^.  ■  9,  as  in  Equation 
(8.2),  the  number  of  variables  are  58320  V3's  ,  9720  a's,  and  230A0  Y's,  giving 
a  total  of  91,080  variables.  If  we  count  9720  V2's  and  11,520  V3's  instead  of 
V3's,  the  total  number  of  variables  is  54,000.  Using  the  CRAY  computer,  150 
major  Iterations  were  done  In  50  minutes  CPU  time. 
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When  we  used  the  CRAY  computer  at  AFWL,  the  case  of  nc  =  9  was  done,  but 
the  nc  =  15  case  was  not,  due  to  the  lack  of  memory  space.  The  VAX/VMS 
computer  has  a  virtual  memory  system  and  could  calculate  the  nc  =  15  case. 

f.  Profile.  The  density  profile  of  the  boundary  is  shown  in  Figure  6. 
The  <  numbers  are  those  in  Figure  2.  The  k  *  2  point  is  the  occupied  lattice 
point  in  the  homogeneous  end-phases.  The  <  =  6  and  7  points  are  the  occupied 
points  on  the  left  end  and  the  unoccupied  points  on  the  right  end. 
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Figure  6.  Density  of  atoms  on  each  sublattice,  at 
kT/e  *  A  and  w/e  ■  -4.  The  density  of 
the  homogeneous  phase  is  0.843. 
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Conclusion 


A  model  of  a  three-dimensional  grain  boundary  has  been  calculated  using 
statistical  mechanics  for  the  first  time. 

The  cluster  variation  method  (CVM)  is  capable  of  treating  such  a  three- 
dimensional  problem.  A  set  of  simultaneous  equations  of  a  large  number  of 
variables  (54,000  variables  or  91,080  variables,  depending  on  how  they  are 
counted)  can  be  solved  using  the  natural  iteration  method  (NIM). 

The  CRAY  computer  can  calculate  the  boundary  about  ten  or  more  times  faster 
in  CPU  time  than  the  VAX/VMS  computer. 
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Abstract  -  The  characteristics  of  equilibrium  first-order  phase  changes  of  the 
grain  boundaries  themselves  are  such  that  It  leads  to  a  definition  of  grain  bound¬ 
ary  phases  In  which  smoothly  curving  boundaries  are  of  the  same  phase.  Different 
grain  boundary  phases  coexist  at  facet  edges  and  corners.  The  phase  rule,  phase 
diagrams,  and  some  phase  change  mechanisms  are  developed.  For  a  wide  variety  of 
problems  orientation  of  the  normal  is  shown  to  be  analogous  to  composition  in 
ordinary  three-component  systems.  The  role  of  symmetry  in  modifying  the  phase 
rule  and  in  sectioning  phase  diagrams  is  explored.  Reports  of  boundary  phase 
changes  are  re-examined  critically. 

Introduction  -  Transformations  of  grain  boundary  structures  have  been  Invoked 
to  explain  observations  of  discontinuities  of  grain  boundary  properties  as  temper¬ 
ature,  composition  and  even  thermal  history  are  varied  [1-9].  This  explanation 
has  led  to  the  concept  of  grain  boundary  phase  transformations  which  assumes  that 
grain  boundaries  are  two-dimensional  phases,  that  there  can  be  several  distinct 
such  phases  each  having  a  free  energy,  and  that  there  will  be  a  phase  change 
whenever  these  free  energy  curves  cross.  At  such  transitions,  properties  can 
change  dlscontinuously  but,  as  has  been  pointed  out,  there  are  limits  [2,10,11]. 

The  boundary  free  energy  can't  be  discontinuous  as  has  been  reported  [3]. 

Phase  changes  are  of  single  crytal  surfaces  are  well  documented  [12-13].  Careful 
studies  on  planar  surfaces  reveal  their  two  dimensional  space  group  symmetry. 
Discontinuities  in  properties,  such  as  adsorption,  have  not  only  been  documented, 
but  symmetry  changes,  which  necessarily  mark  a  phase  change  that  could  be  higher 
order,  occur  frequently. 

Grain  boundaries  have  rarely  been  studied  with  such  care.  Even  the  best  studies 
have  dealt  with  curved  boundaries  and  since  properties  of  boundaries,  including 
their  symmetry,  depend  on  orientation,  the  question  of  what  is  a  grain  boundary 
phase  needs  to  be  answered  before  we  can  ask  if  there  has  been  a  change  of  phase. 

In  this  paper  we  defer  this  question.  We  assume  that  there  are  indeed  several 
stable  or  metastable  structures  of  a  grain  boundary  for  each  orientation.  We 
compare  these  to  find  the  most  stable  structure  or  combination  of  structures. 

From  this  comes  the  surprising  result  that  only  a  few  specific  orientations  can 
transform  congruently,  l.e.  without  a  change  of  shape.  Almost  all  grain  boundary 
phase  transitions  involves  facetting  to  different  orientations.  As  a  consequence, 
a  natural  definition  of  grain  boundary  phases  and  coexisting  phases  at  equili¬ 
brium  presents  Itself.  The  several  types  of  phase  transitions  then  are 
re-examined.  We  include  wetting  transition  In  which  a  phase  not  stable  In  bulk 
forms  along  a  grain  boundary  but  we  exclude  the  bulk  phase  changes  that  originate 
at  grain  boundary. 

We  begin  by  specifying  the  state  variables  of  grain  boundaries,  examine  conditions 
for  equilibrium  and  then  develop  the  properties  of  diagrams  of  state  or  phase 
diagrams. 

6ra in  Boundary  State  Variables  -  As  with  three-dimensional  phases,  there  are  a  num- 
ber  of  variables  that  define  the  state  of  a  planar  grain  boundary.  The  structure 
and  properties  of  an  equilibrated  boundary  are  fully  specified  when  these  variables 
are  fixed.  Because  the  boundary  Is  free  to  exchange  matter  and  energy  with  the 
grains.  It  acts  as  an  open  system.  As  a  result,  temperature  and  chemical  potential 
are  the  natural  thermodynamic  variables  for  describing  the  reservoir.  Gibbs'  proof 
[14]  that  temperature  and  chemical  potential  are  constant  throughout  a  system  of 
surfaces  and  phases  is  readily  extended  to  Include  grain  boundaries  [2,14,15]. 

Since  grain  boundary  temperature  and  chemical  potentials  are  specified  by  temper¬ 
ature,  pressure,  and  composition  of  a  grain,  these  could  equally  well  serve  as 
state  variables  of  the  grain  boundary.  With  this  latter  set  the  grain  boundary 
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remains  an  open  system  as  long  as  the  grain  size  is  large  enough  that  changes  in 
adsorption  will  not  affect  the  grain  composition.  The  alternate  choice  of  a  compo- 
slitional  state  variable,  the  level  of  adsorption,  has  many  disadvantages.  It  is 
not  as  easily  measured  and  in  an  equilibrated  polycrystalline  specimen  varies  from 
place  to  place  on  the  boundaries. 

In  addition  to  these  standard  thermodynamic  variables,  the  state  of  the  boundary 
depends  on  geometric  variables.  Three  angular  variables  specify  R,  the  misorien- 
tation  between  the  two  abutting  grains.  Two  angular  variables  specify  N,  the 
orientation  of  the  grain  boundary  normal.  There  are  four  more  translational  varia¬ 
bles  [16]  that  we  will  ignore  because  they  are  usually,  but  not  always,  allowed  to 
relax  to  equilibrium.  In  some  experiments  the  angular  variables  are  also  allowed 
to  relax  to  equilibrium  [17],  Grains  can  rotate  and  boundaries  can  facet  to  other 
orientations.  We  will  specifically  assume  that  R  and  the  average  of  N  are 
specified. 

The  total  number  of  grain  boundary  state  variables  is  therefore  six  plus  the 
number  of  chemical  components.  There  are  five  angular  variables  (R  and  N),  temper¬ 
ature,  and  either  one  chemical  potential  for  each  of  the  components  or  pressure  and 
the  composition.  Area  is  the  extensive  variable.  We  will  mostly  specify  area 
projected  along  N.  This  is  the  area  that  is  fixed  if  we  specify  a  planar  curve  to 
be  the  perimeter  of  the  boundary.  The  boundary  then  is  able  to  find  a  minimum  free 
energy  at  the  fixed  projected  area  even  if  that  minimum  includes  facetting  or 
zigzagging  to  orientations  of  lower  energy  [18], 

Grain  Boundary  Equilibrium  -  The  appropriate  free  energy  is  one  that  is  minimized 
at  equilibrium  under  the  conditions  specified  by  the  method  of  the  observation. 

We  fix  the  perimeter  of  the  grain  boundary  to  be  a  plane  curve  and  permit  matter 
and  energy  to  be  exchanged  with  the  abutting  grains  of  fixed  temperature  and  either 
fix  chemical  potentials  or  composition  and  pressure.  For  grain  boundary,  as  well 
as  surface  equilibrium,  it  is  the  excess  in  E-TS-IuiNj  per  unit  area  projected 
along  N  that  is  to  be  a  minimum  [15,18,19].  This  excess  free  energy  per  unit 
area  has  been  variously  described  by  the  symbols  o  or  y.  We  use  y. 

Equilibration  of  a  grain  boundary  involves  reducing  y  to  its  lowest  value. 

This  may  involve  exchanges  of  energy  and  matter  with  the  grains,  structural  rear¬ 
rangement,  and  facetting.  For  the  purposes  of  this  paper  we  consider  not  just  the 
lowest  value  of  y  but  several  grain  boundary  structures  that  are  in  metastable  or 
unstable  equilibrium.  These  other  boundaries  are  at  stationary  values  of  y  that 
are  local  minima. 

Such  other  structures  are  a  familiar  concept  in  metastability.  They  provide  a 
vehicle  for  the  description  of  first  order  phase  changes  in  which  quite  different 
structures  become  the  stable  phase  as  the  state  variables  are  varied  and  give 
rise  to  a  phase  diagram  whose  dimension  equals  the  number  of  state  variables, 
which  maps  the  bounds  to  the  stability  of  a  grain  boundary  phase. 

First  Order  Grain  Boundary  Phase  Transformations  -  The  simplest  way  to  examine 
grain  boundary  phase  transformations  is  to  image  two  or  more  different  grain  bound¬ 
ary  structures  and  to  examine  graphically  how  their  y  values  vary  with  one  of  the 
nongeometric  state  variables.  Thermodynamics  requires  that  for  these  variables 
each  of  the  y  plots  be  continuous.  Let  us  examine  a  range  of  temperature  over 
which  two  such  hypothetical  y  curves  cross.  For  two  forms  of  a  fluid  surface  in  a 
one  component  constant  pressure  system,  there  would  be  a  phase  change  between  these 
two  forms  exactly  where  the  two  curves  cross.  This  situation  is  identical  to  a 
first-order  one-component  phase  change  in  three  dimensions  with  Gibbs  energy 
Instead  of  y  as  the  ordinate. 

To  demonstrate  that  this  simple  example  Is  misleading  for  solid  boundaries,  we 
turn  to  a  familiar  problem  in  another  field.  It  is  well  known  that  phase  changes 
In  a  system  with  fixed  average  composition  do  not  usually  occur  simply  at  the  place 
where  the  molar  Gibbs  free  energies  of  the  phases  cross  [20-22].  Figure  1  shows  a 
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pair  of  Intersecting  isothermal  constant  pressure  free  energy  curves  in  a  binary 
system.  The  common  tangent  represents  a  two-phase  equilibrium  of  coexisting 
phases,  differing  in  composition,  at  a  lower  free  energy  than  either  curve  in  that 
composition  range.  The  loci  of  intersections  of  free  energy  curves  in  such  a 
system  do  not  mark  phase  changes.  Because  only  the  average  composition  is  spec¬ 
ified,  the  system  is  free  to  seek  equilibrium  with  a  lower  free  energy  in  a 
combination  of  phases  of  differing  compositions.  This  is  the  essential  distinction 
between  the  constructions  in  the  one  and  two  component  systems:  both  phases  must 
have  the  same  temperature,  while  only  the  average  composition  is  specified. 

Equilibrium  is  expressed  by  a  convexity  theorem,  which  requires  that  the  isothermal 
molar  free  energy  vs.  composition  curve  be  convex  as  seen  from  below.  This  convex 
curve,  consisting  of  segments  of  curves  and  tangent  lines,  is  depicted  in  Figure  1. 
The  curved  portions  represent  single  phases,  the  tangent  lines  represent  the  free 
energies  of  mixtures  of  phases  whose  compositions  are  those  of  the  contacts.  This 
construction  eliminates  not  only  any  concave  bulges  but  also  the  intersections  of 
free  energy  curves.  As  temperature  and  pressure  change  the  free  energy  curves,  the 
compositions  of  the  two-phase  region  shifts.  These  shifts  are  graphically  repre¬ 
sented  on  phase  diagrams.  Of  special  interest  are  conditions  where  two  free 
energycomposition  curves  touch  without  crossing.  These  represent  congruent  phase 
changes,  where  two  phases  of  the  same  composition  are  in  equilibrium. 

In  multicomponent  systems  the  free  energy-composition  plot  at  constant  temperature 
and  pressure  is  a  hypersurface.  Equilibrium  requires  that  this  hypersurface  also 
be  convex  and  that  concave  bulges  be  eliminated  by  hyperplanes.  The  phase  rule 
[14,20,21]  is  an  expression  of  how  many  phases  can  coexist  and  has  a  very  simple 
geometrical  interpretation  in  terms  of  these  tangents.  It  must  be  amended  when 
symmetry  is  present,  as  in  mixtures  of  isotopes  or  enantiomorphs  [23].  Symmetry 
of  the  free  energy  hypersurface  can  permit  more  simultaneous  tangent  contacts. 

An  important  aspect  of  phase  diagrams  concerns  which  of  two  kinds  of  variables  are 
chosen  as  axes  [24],  Quantities  such  as  composition  need  not  be  the  same  in 
coexisting  phases.  Hence,  they  lead  to  gaps  in  compositions  that  are  multiphase 
regions  and  to  tie  lines,  tie  triangles,  etc.,  across  these  gaps  to  indicate  the 
composition  of  coexisting  phases.  Quantities  such  as  temperature,  pressure,  and 
chemical  potentials  have  the  same  value  in  all  coexisting  phases.  If  they  are 
axes  on  a  phase  diagram,  the  tie  lines  must  be  perpendicular  to  these  axes.  If 
all  axes  are  of  this  type  no  gaps  can  occur.  Multidimensional  phase  diagrams 
sectioned  perpendicular  to  such  axes  are  easier  to  interpret  because  coexisting 
phases  always  lie  within  such  a  section. 

We  now  imagine  two  different  grain  boundary  phases  that  are  represented  by  two 
Intersecting  y  curves.  For  the  present  each  curve  is  assumed  to  be  a  continuous 
function  of  all  the  state  variables  including  N.  Does  the  locus  of  points  of 
Intersection  mark  the  condition  of  phase  transition?  Because  of  the  open  system 
conditions,  the  temperature  and  chemical  potentials  of  all  coexisting  grain 
boundary  phases  are  the  same.  Thus  we  must  compare  grain  boundary  phases  at  the 
same  temperature  and  chemical  potentials.  This  is  true  even  if  we  use  grain 
compositions  as  our  variables.  The  geometric  variables  require  special  examina¬ 
tion.  We  may  try  to  Impose  a  misorientation  R,  but  grain  boundaries  are  known 
[25]  to  dissociate  Into  several  grain  boundaries  whose  combined  misorlentations 
are  equivalent  to  R.  The  orientation  N  has  the  character  that  only  its  average  is 
usually  specified  but  the  boundary  Is  free  to  facet  or  zigzag. 

There  are  several  constructions  that  are  useful  for  examining  the  N  dependence* 
fterrlng  [18]  derived  a  tangent  sphere  construction  for  examining  stability 
which  Frank  [26]  simplified  by  inverting  y  in  which  the  test  sphere  became  a 
tangent  plane.  This  was  then  used  by  Hoffman  and  Cahn  [19,27]  to  construct  £  a 
vector  function  of  y.  Both  of  the  last  two  lend  themselves  readily  to  examining 
the  N  dependence. 

Frank  took  the  reciprocal  of  y  and  showed  that  a  plot  of  y"*  In  polar  coordinates 
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had  to  be  convex;  i.e.,  any  tangent  plane  to  yl  had  to  stay  outside  the  yl 
plot.  Inward  bulges  were  to  be  filled  in  with  tangent  planes  representing  hill  an< 
valley  structures  (Fig.  1).  The  geometrical  convexity  condition  for  this  three- 
dimensional  surface  is  mathematically  identical  to  the  condition  encountered  in 
three-component  phase  equilibria.  The  behavior  of  N  in  surface  and  grain  boundary 
phase  transitions  is  identical  to  the  behavior  of  two  compositional  variables  in  a 
three-component  system. 


Fig.  1  :  The  same  tangent  construction  Is  used  for  the  coexistence 
of  phases,  using  a  cartesian  plot  of  molar  free  energy  vs.  composi¬ 
tion,  and  of  facets,  using  a  polar  plot  of  y-^N).  The  tangent 
eliminates  all  concave  portions,  whether  caused  by  a  bulge  in  a 
single  curve  or  by  intersecting  convex  curves. 

The  development  of  the  £  vector  leads  to  identical  conclusions.  The  £  vector  is 
constructed  from  y  [19,27].  The  loci  of  £  vectors  forms  a  surface  whose  inner 
envelope  is  a  closed  convex  surface  whose  shape  is  that  of  a  grain  of  fixed  volume 
embedded  within  another,  misoriented  by  R,  with  minimum  total  grain  boundary  free 
energy.  The  £-plot  represents  the  actual  orientations  present.  Orientations 
that  are  unstable  with  respect  to  facetting  are  absent,  and  graphically  seen  as 
corners  and  edges.  The  edges  of  any  £  plot  mark  two  coexisting  orientations; 
corners  mark  coexistence  of  more  than  two  orientations.  If  there  were  two  kinds 
of  grain  boundary  structures  with  intersecting  y,  there  would  be  two  intersecting 
£  plots  (Fig.  2),  whose  inner  envelope  would  necessarily  have  a  sharp  edge  indica¬ 
ting  that  certain  orientations  would  always  be  absent  at  equilibrium.  The  coexist¬ 
ing  phases  would  meet  at  an  edge  in  the  £-plot.  They  would  display  an  edge  with 
the  same  orientations  in  a  real  specimen  at  equilibrium. 

Convex  portions  of  £-plots  outside  the  inner  envelope  are  metastable,  and  facet 
by  nucleation  of  patches  of  a  new  orientation.  Portions  that  are  not  convex  are 
unstable  with  respect  to  continuous  reorientation.  Concave  bulges  In  Figure  1 
lead  to  self  intersections  in  Figure  2  with  convex  and  concave  portions  separated 
by  spinodals  [27-29], 

Except  for  congruent  phase  changes,  and  the  unlikely  possibility  that  y  is  not 
A  function  of  N,  all  equilibrium  phase  changes  involve  the  formation,  elimination, 
or  other  changes  in  edges  and  corners  of  a  grain  boundary.  Equilibria  between 
these  different  structures  would  almost  always  occur  as  observed  edges  and  corners 
In  a  hill  and  valley  structure.  Facetting  can  also  occur  if  we  assume  a  single  y 
curve  with  a  concave  bulge.  Such  free  energy  bulges  are  analogous  to  those  encoun¬ 
tered  In  miscibility  gap  systems  that  have  critical  points.  The  coexisting  phases 
are  different  phases.  At  such  critical  points  the  phases  become  Identical.  The 
analogy  holds  for  grain  boundary  phases.  75 
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We  now  propose  to  define  grain  boundary  phases  such  that  all  equilibrium  edges  and 
corners  represent  coexisting  phases.  This  Includes  all  coexistence  based  on  the 
assumption  of  several  metastable  structures  with  the  same  N.  In  this  definition 
curved  grain  boundary  without  singularities  In  Its  shape  Is  all  of  the  same  phase. 
The  structure  of  this  grain  boundary  phase  varies  as  the  orientation  changes,  just 
as  the  structure  of  a  chemical  phase  varies  over  Its  composition  range.  Phases 
differing  by  a  first-order  phase  change  coexist  at  edges  and  corners  (except 
when  It  Is  a  congruent  phase  change).  Singularities  in  shape  other  than  dlscontln 
ultles  In  orientation  represent  higher  order  phase  changes. 


Fig.  2  :  The  Inner  envelope  of  a 
£  plot  shows  the  actual  ori¬ 
entations  of  stable  boundaries. 
Orientation  discontinuities  (edges) 
represent  coexisting  orientations 
(facets).  Intersecting  £  plots 
necessarily  lead  to  orientation 
discontinuities.  Self  Intersecting 
portions  of  the  £  plot  result 
from  bulges  in  y"1  plots,  and 
show  "ears"  with  Inelastable 
convex  portions  and  unstable  con¬ 
cave  portions  separated  by  spl nodal s 


Sometimes  the  coexisting  grain  boundary  phases  have  Identical  properties  by  symme¬ 
try.  This,  too,  has  an  analogy  In  three-dimensional  phase  diagrams.  Systems  with 
enantlomorphs  have  different  phases  that  are  mirror  images  of  each  other  with 
Identical  properties.  As  In  the  case  of  enantlomorphs,  symmetry  leads  to  modifica¬ 
tions  of  the  phase  rule  for  grain  boundary  phases  [23]. 

We  began  by  assuming  the  existence  of  two  different  structures  with  two  dif¬ 
ferent  y  curves  and  concluded  that,  except  for  the  congruent  orientation  where  two 
different  planar  boundaries  transform  Into  each  other,  most  phase  changes  Involve 
the  formation  of  patches  of  grain  boundary  with  different  orientations. 

We  also  began  by  assuming  that  y  was  a  continuous  function  of  N.  This,  too,  is  not 
essential.  Frank's  y’1  and  our  £  plot  are  continuous  convex  figures  even  If  y  Is 
Is  not.  Planes  or  ruled  surfaces  alternate  with  doubly  curved  portions.  In  the 
polyhedra  as  the  curved  portion  have  shrunk  to  points  on  the  y"*  plot  and  have 
become  the  planes  on  the  £  plot. 

The  convex  y“*  plots  are  continuous  functions  of  N.  The  Gibbs  adsorption  equa¬ 
tion  Insures  that  y  is  a  continuous  function  of  R  and  the  chemical  potentials. 
Stability  consideration  requires  that  y  either  be  continuous  In  R  or  dissociate 
Into  boundaries  whose  combined  y  Is  continuous  In  R.  Thus  y  Is  continuous  through 
phase  transitions  whether  the  transition  Is  congruent  or  proceeds  by  facetting,  or 
which  of  the  state  variables  Is  varied.  Reports  that  y  be  discontinuous  Is  con-  • 
trary  to  thermodynaml c  stability  principles.  It  cannot  be  explained  by  any  arti¬ 
fact  such  as  impurities  [3,5]  which,  after  all,  are  just  chemical  components. 
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Many  correlations  follow  from  this  rigorous  mathematical  analogy  between  orienta¬ 
tion  of  surfaces  and  composition  in  a  three  component  system.  These  will  be 
developed  in  the  next  section  as  part  of  an  examination  of  all  the  state  variables 
and  the  full  phase  diagram. 

Grain  Boundary  Phase  Diagrams  -  We  next  examine  a  space  of  all  6+C  state  variables, 
where  C  is  the  number  of  chemical  components.  Five  axes  are  angular  variables, 
three  to  specify  R  and  two  to  specify  N.  These  specify  temperature  the  chemical 
potentials  or  activities  of  the  components.  Each  point  in  this  hyperspace  repre¬ 
sents  a  state  of  a  grain  boundary.  Some  regions  represent  stable  planar  boundaries 
of  a  particular  phase,  other  regions  represent  the  fact  that  a  boundary  with  these 
state  variables  is  decomposed  into  segments  of  several  boundaries,  each  segment 
being  of  a  different  phase. 

The  examination  of  this  hyperspace  Is  greatly  simplified  by  the  realization  that 
any  set  of  coexisting  grain  boundary  phases  share  the  same  temperature  and  chem¬ 
ical  potentials.  Hence  grain  boundary  phase  diagrams  can  be  sectioned  at  constant 
temperature  and  chemical  potentials  with  the  knowledge  that  all  tie  lines  (or  tie 
figures  of  higher  connectivity),  connecting  points  that  represent  the  coexisting 
boundary  phases,  lie  entirely  within  that  section.  This  section  still  has  the  five 
dimensions  of  the  R  and  N  axes,  but  the  6+C  dimensional  diagram  can  now  be  thought 
of  as  an  evolution  of  five-dimensional  diagrams  as  T  and  the  chemical  potentials 
are  varied. 

The  five-dimensional  diagram  can  be  sectioned  further.  If  all  facets  have  the  same 
misorientation  R  across  them,  then  all  tie  figures  lie  within  a  constant  R  section. 
This  is  a  two-dimensional  section  with  N  as  its  axes.  Such  N  diagrams  have  the 
features  of  an  isothermal  constant-pressure  section  of  a  ternary  phase  diagram 
plotted  onto  the  surface  of  a  unit  sphere  instead  of  an  equilateral  triangle. 

Although  R  can  be  Imposed  by  creating  a  bl crystal  with  that  misorientation,  the 
system  will  dissociate  this  boundary  into  several  with  a  combined  R  equal  to  the 
original  misorientation  If  that  will  lower  the  energy.  One  such  dissociation  has 
been  demonstrated  [25].  When  a  high  energy  boundary  has  the  combined  misorien¬ 
tation  of  two  known  low  energy  boundaries,  such  as  twins  or  special  coincidence 
site  boundaries,  such  boundaries  dissociate.  The  resultant  coexistent  boundaries 
have  R  values  different  from  the  imposed  R  and  therefore  do  not  lie  In  a  constant 
R  section.  The  twin  part  of  the  dissociated  boundary  has  been  observed  to  facet 
[25].  In  this  case  the  tie  lines  cannot  be  counted  on  to  be  perpendicular  to  any 
axis.  Unless  there  are  some  special  symmetries,  the  representation  of  the  phase 
coexistence  requires  the  full  five-dimensional  space. 

The  construction  rules  for  coexistence  of  dissociated  boundaries  differs  from  that 
of  the  usual  phase  diagram.  The  combined  R  is  conserved  in  a  dissociation,  where¬ 
as  the  average  N  or  composition  is  conserved  in  facetting  or  phase  separation.  As 
a  result,  the  limits  of  regions  that  represent  undissociated  boundaries  do  not 
represent  R  values  of  the  components  of  a  dissociated  boundary.  No  method  of  re¬ 
presenting  dissociations  on  the  diagram  has  been  suggested. 

The  sections  of  the  C+6  dimensional  phase  diagram  parallel  to  temperature  or 
chemical  potential  axes  may  have  tie  lines  at  an  angle  to  the  section.  This  is  a 
familiar  complication  in  "vertical"  sections  of  ternary  phase  diagrams  in  which 
temperaturfe  or  pressure  are  the  vertical  axis.  If  there  Is  no  dissociation,  such 
sections  should  look  quite  similar  to  such  vertical  sections. 

Vertical  sections  containing  all  tie  lines  exist  for  ternaries.  These  sections 
are  true  binaries.  They  are  sometimes  called  quasibinaries.  Symmetry  is  an  im¬ 
portant  guide  for  finding  such  vertical  sections  for  grain  boundaries.  (Vertical 
sections  of  surface  phase  diagrams  have  been  given  [12,13]  for  symmetric  orien¬ 
tations.)  77 


The  Phase  Rule  -  Many  years  ago,  Oefay  [31,32]  proposed  a  surface  and  Interface 
phase  rule  for  systems  with  several  bulk  fluid  phases  and  allowed  for  possibility 
of  a  multitude  of  interfaces  between  pairs  of  different  phases.  He  dealt  only 
with  the  usual  C+2  state  variables.  His  phase  rule  should  not  be  applied  to 
-solids  because  It  does  not  take  into  account  the  multitude  of  solid  phenomena, 
e.g. ,  orientation  effects.  What  we  develop  here  Is  a  localized  grain  boundary 
phase  rule  for  a  single  phase  polycrystal  line  solid.  It  deals  with  how  many  grain 
boundary  phases  can  coexist  at  a  point  on  a  grain  boundary  and  whether  this  coexis¬ 
tence  remains  If  we  change  the  state  variables.  The  rule  Is  easily  modified  If 
there  are  several  bulk  phases  but  it  remains  a  local  phase  rule.  The  phase  rule 
Is  a  useful  concept  for  phase  diagrams  of  first  order  transitions.  The  derivation 
of  a  phase  rule  assumes  that  all  grain  boundary  phases  have  different  fundamental 
equations  for  how  y  depends  on  the  state  variables.  This  Is  frequently  not  the 
case.  When  high  symmetry  grain  boundaries  facet,  they  will  facet  Into  symmetry 
related  orientations,  and  the  phase  rule  must  be  amended  for  symmetry.  Fortun¬ 
ately  the  y"*  and  £  constructions  are  geometrical  constructions  completely 
consistent  with  the  phase  rule  and,  as  will  be  shown,  these  permit  us  to  specify 
what  happens  at  the  exceptions. 

When  all  grain  boundary  phases  are  different,  the  phase  rule  specifies  that  the 
degrees  of  freedom  F  is  given  by 

F  *  C  -  P  +  7  (1) 

where  C  Is  the  number  of  chemical  components,  and  P  the  number  of  grain  boundary 
phases  all  nonequivalent.  When  there  is  more  than  one  bulk  phase  the  phase  rule 
becomes 

F  *  C  -  (P  +  PB)  +  8  (2) 

where  PB  is  the  number  of  bulk  phases. 

If  R,  T,  and  the  chemical  potentials  are  specified  (or  equivalently  pressure  and 
the  composition),  we  obtain  the  rule  for  an  N  diagram,  again  for  nonequivalent 
phases. 

F  *  3  -  P  (3) 

r 

which  Is  the  same  as  for  a  three-component  phase  diagram  at  constant  T  and 
pressure.  We  examine  these  N  diagrams  with  special  emphasis  on  the  symmetry 
effects.  Because  crystal  surfaces  share  all  the  same  state  variables  with  grain 
boundaries  except  R,  the  properties  of  this  section  apply  equally  to  surfaces. 

We  have  given  a  more  complete  treatment  for  surfaces. 

The  N  Diagram  -  In  order  to  examine  the  properties  of  an  N  diagram,  we  assume 
that  all  tie  lines  lie  within  that  section.  R  Is  assumed  fixed  and  stable  with 
respect  to  dissociation.  We  then  know  a  great  deal  about  the  properties  of 
such  a  section  for  It  Is  equivalent  to  a  Isothermal  constant  pressure  section  of 
a  three-component  phase  diagram  plotted  on  the  surface  of  the  unit  sphere  instead 
of  an  equivalent  composition  triangle. 

We  first  examine  phase  coexistence  when  all  phases  are  nonequivalent.  In  an  N 
diagram,  because  R,  T,  and  the  chemical  potentials  are  specified,  equation  (3) 
applies.  Single-phase  regions  have  the  two  degrees  of  freedom  that  permit  N  to 
range.  Smoothly  curved  boundaries  In  an  experiment  plots  are  single  phases.  So 
are  smoothly  curved  regions  of  a  £  plot,  as  are  those  curved  convex  regions  on 
a  y"1  plot.  The  single  phase  on  the  N  diagram  can  occupy  the  entire  unit 
sphere,  or  It  can  be  bounded  by  phase  changes.  Two-phase  coexistence  Is  repre¬ 
sented  by  two  curves  (F  <*  1)  connected  by  tie  lines  representing  coexisting 
orientations.  The  space  on  the  diagram  between  the  curves  are  forbidden  orien¬ 
tations,  edges  In  an  experiment  and  on  a  £-p1ot.  The  tie  lines  and  coexistent 
phases  are  formed  by  a  rolling  tangent  plane  on  the  y“*  plot.  Tie  lines  are  75 


great  circle  segments,  whose  poles  are  directions  of  the  tangents  to  the  edges 
they  represent. 

These  two-phase  regions  may  form  smooth  closed  curves.  They  may  also  terminate 
at  critical  point  where  the  two  phases  have  become  identical  or  at  three  phase 
coexistence.  On  a  5-plot  critical  points  are  end  points  of  edges  In  which 
the  crease  terminates  and  Is  smooth  thereafter.  It  Is  represented  by  a  termina¬ 
ting  orientation  gap  on  an  N  diagram.  Three-phase  coexistence  is  a  corner  in  the 
experiment  or  on  a  5-plot,  a  tangent  contacting  the  y"1  plot  at  three 
points.  It  is  represented  by  a  tie  triangle  on  the  N-diagram. 

To  obtain  congruent  points  or  “more  than  three-phase"  coexistence,  one  or  more  of 
the  other  state  variables  must  be  changed,  or  symmetry  has  to  be  present. 

Both  y  as  a  function  of  N  and  the  N-diagram  must  obey  the  Wulff  group  point 
symmetry  W  Imposed  by  R  [33].  For  surfaces  this  point  symmetry  Is  that  of  the 
crystal.  For  grain  boundaries,  it  is  derived  [16,34,35]  from  the  point 
symmetry  of  the  grains  and  R.  The  set  W  of  resulting  point  groups  include  the 
usual  32  point  groups  plus  14  additional  with  8  and  12  fold  axes,  a  total  of  46. 

We  call  these  latter  point  groups  octagonal  and  dodecanogal. 

The  effect  of  symmetry  Is  most  easily  seen  using  the  5-plots.  Nothing  Is 
changed  for  single  and  two-phase  coexistence,  which  are  still  represented  by 
surfaces  and  curved  edges.  Only  two  surfaces  can  meet  at  an  edge  of  a  convex  Inner 
envelope.  More  than  three  surfaces  can  meet  at  corners.  We  may  use  crystal¬ 
lographic  form  theory  [36]  to  specify  how  many. 

For  the  32  classical  point  groups,  there  are  47  forms  resulting  from  Intersections 
of  planes  that  are  equivalent  by  symmetry.  Many  are  prisms  and  pinacolds  composed 
of  planes  parallel  to  an  axis,  that  do  hot  form  corners.  The  remainder  have 
corners  where  three  or  more  equivalent  planes  meet.  Corners  also  result  from 
combined  forms  In  which  planes  that  are  not  all  equivalent  intersect  at  a  point. 
There  are  additional  forms  that  originate  In  the  octagonal  and  dodecagonal  point 
groups. 

There  are  two  classes  of  forms.  There  Is  one  general  form,  in  which  the  planes 
haye  general  indexes,  for  each  of  the  46  point  groups  except  1.  There  are  In 
addition  special  forms  for  special  Miller  Indexes. 

We  examine  all  the  forms  for  corners.  Table  I  gives  the  number  of  facets  along 
the  axis  for  tetragonal,  hexagonal,  octogonal  and  dodecagonal  forms  for  general 
Indexes  and  the  two  kinds  of  special  Indexes  that  do  not  give  pinacolds  or  prisms. 
This  number  Is  equal  to  the  number  of  coexisting  facets  if  the  basal  plane  facets. 
Figure  3  shows  the  dlhexagonal  dipyramid  and  the  trigonal  dlpyramid.  These  are 
general  forms  In  the  hexagonal  (r  =  6)  system  for  6/mmm  and  b  respectively  show 
12  and  3  sided  corners.  If  the  basal  plane  were  to  facet  Into  planes  with  general 
Indexes,  there  would  be  12  and  3  facets  if  the  Wulff  group  were  6/mmm  or  C  respec¬ 
tively.  It  Is  Interesting  to  note  that  the  number  of  faces  either  remains  the  same 
or  decreases  if  faceting  occurs  at  a  special  orientation.  The  Indexes  of  the 
facets,  of  course,  refer  to  those  of  the  dichromatic  pattern  or  the  two  interpene¬ 
trated  crystals  and  not  to  those  of  either  crystal. 

Table  II  gives  the  number  of  facets  In  the  cubic  system  If  either  the  cube  plane  or 
the  octahedral  plane  facets.  The  numbers  divisible  for  3  are  those  for  the  octahe¬ 
dral  plane.  Whenever  there  are  parenthesis  around  them  (111)  and  (111)  give  dif¬ 
ferent  answers.  In  contrast  to  Table  I,  the  number  of  facets  Is  often  larger  for 
special  orientations.  All  faceting  of  basal  planes  in  the  trigonal  system  yields 
3,  except  for  hkll  and  hh2hl  In  point  groups  32  and  3m  which  yield  6.  Four-sided 
corners  occur  along  the  2-fold  axis  in  mm2  and  along  all  three  axes  in  mmm  for 
special  indexes  only.  All  other  four-sided  corners  arise  from  mirror  planes  or 
combined  forms.  79 
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Fig.  3  :  Forms  of  5  plots  must  conform  to  the  symmetry  of  the  Wulf£ 
group.  The  basal  plane  will  facet  into  3  and  12  orientations  for  6 
and  6/mmm  respectively.  Note  also  the  meeting  of  four  orientations 
at  the  mirror  plane  in  6/mmm. 
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Table  1. 

The  number  of  equivalent  facets  replacing  the  basal 

plane,  listed  by  point 

group 

and  facet  indexes  when 

the  rotation  axes  are 

r-fold  with  r  *  4,  6,  8, 

and  12 

Point  Group 

Facet  Index 

hkl 

hhl 

hOl 

* 

hkll 

hohl 

hh2hl 

r 

r 

r 

r 

r 

- 

r/2 

r/2 

r/2 

r/m 

r 

r 

r 

r22 

r 

r 

r 

rmm 

2r 

r 

r 

rm2 

r 

r/2 

r 

r/mmm 

2r 

r 

r 

Table  II.  The  number  of  equivalent  facets  replacing  the  cube  and  octahedral 
planes  for  cubic  systems. 
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Four-sided  corners  arise  along  mirror  planes  whenever  an  edge  Intersects  a  mirror 
plane.  This  can  be  seen  along  the  basal  plane  of  the  dlhexagonal  dipyramid.  It 
also  occurs  In  combined  forms  In  systems  with  symmetry  as  low  as  m  (Fig.  4). 

Although  these  multiphase  equilibria  range  from  P  ■  3  to  P  »  24  they  occur  In  N 
diagrams  without  requiring  more  degrees  of  freedom  than  for  P  ■  3.  If  a  tangent 
plane  contacts  the  y-1  surface  at  three  equivalent  points.  It  contacts  all 
points  that  are  equivalent  by  symmetry.  If  a  tangent  plane  contacts  two  points  on 
one  side  of  a  mirror,  it  can  contact  either  none  or  two  points  on  the  other  side. 

If  we  use  the  other  degrees  of  freedom  to  adjust  the  relative  y  values  of  different 
phases,  we  can  achieve  multiple  contacts  in  combined  forms.  Each  additional  form 
requires  a  degree  of  freedom,  just  as  It  does  when  all  phases  are  different  in 
equation  1.  No  more  than  three  facet  orientations  are  ever  needed  to  reach  a 
minimum  free  energy  [37]  for  an  undissociated  boundary,  or  for  each  component 
boundary  In  a  dissociated  one.  When  the  phase  rule  permits  more  facets,  many 
degenerate  structures  are  possible. 

For  the  purpose  of  the  phase  rule,  two  phases  related  by  symmetry  count  as  two 
phases,  while  three  or  more  count  as  three  phases.  M  pairs  of  phases  related  by 
a  mirror  plane  count  as  M+1  phases. 

The  Search  for  Congruent  Transformations  -  Congruent  phase  changes  do  not.  In 
general,  occur In  a  constant  R,  T  and  chemical  potential  section.  This  can  be 
easily  seen  from  the  phase  rule  and  from  the  y**  or  £  plots.  Congruent  points 
occur  on  these  plots  only  when  the  y”1  or  £  plots  of  one  phase  contacts  the 
other  at  a  point  without  Intersecting.  We  therefore  need  to  adjust  one  of  the 
other  state  variables  to  obtain  precisely  this  condition.  The  search  for  congruent 
transformations  has  all  the  aspects  of  a  search  for  a  needle  In  a  6+C  dimensional 
haystack.  Nonetheless  one  congruent  phase  transformation  has  been  found  experimen¬ 
tally  and  theoretically,  making  use  of  known  symmetries.  Congruent  phase  changes 
are  likely  to  occur  at  high  symmetry  orientations  because  y  has  an  extremum  there 
and  the  tangent  contact  requires  the  same  slope  for  both  phase.  But  there  seems  to 
be  no  Increase  In  the  degrees  of  freedom  because  of  symmetry  for  this  transition. 

Gllcksman  and  Void  [-38]  found  a  congruent  grain  boundary  melting  transition  In 
bismuth,  as  the  varied  R.  The  rotation  axis  was  the  trigonal  3  axis.  For  a 
general  rotation  the  dichromatic  pattern  symmetry  Is  3m1  with  the  colored  mirror 
planes  bisecting  directions  In  the  two  crystals  related  by  the  rotation.  The 
Wulff  group  symmetry  for  y  Is  thus  3m.  The  symmetric  tilt  boundaries  studied  lie 
along  these  mirror  planes.  Mirror  planes  by  themselves  do  not  guarantee  a  symmetry 
dictated  extremum,  but  In  a  3m  the  mirrors  are  part  of  a  2/m  axis  that  Is  suffi¬ 
cient.  These  boundaries  are  likely  candidates  for  congruent  transformations  as  R 
Is  varied.  The  nonsymmetric  boundaries  that  do  not  lie  along  the  mirror  planes 
facet  on  cooling  through  the  melting  transition  [39]. 

A  congruent  melting  transition  was  also  found  theoretically  [40]  by  varying  T 
for  a  symmetric  I  *  5  boundary  In  a  two-dimensional  square  lattice  gas  crystal. 
Again  the  required  symmetry  was  present. 

Neither  of  these  studies  provided  clean  examples  of  how  y  should  vary  through  a 
first  order  congruent  transition  l.e.  that  y  be  continuous  but  with  a  discon¬ 
tinuity  In  slope.  In  the  bismuth  experiments,  the  grain  boundaries  were  equi¬ 
librated  with  solid-melt  Interfaces  whose- y  was  measured  and  shown  to  be  almost 
Isotropic.-  The  dihedral  angle  of  the  grain  boundary  groove  changed  dlscontlnuously 
at  the  melting  transition. 

One-dimensional  features  do  not  show  clean  phase  transitions,  only  regions  of  rapid 
change.  The  grain  boundary  between  two  two-dimensional  crystals  is  such  a  one 
dimensional  feature,  and  the  calculated  phase  transition  had  the  expected  gradual 
nature  [40]. 
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The  Controversy  About  Dihedral  Angle  Discontinuities  -  In  1970  Glelter  [3] 
examined  a  trijunction  In  lead  and  reported  a  discontinuity  In  the  dihedral  angles 
as  temperature  was  changed.  Gllcksman  and  Void  [38]  reported  a  discontinuity  In 
the  grain  boundary  melt  groove  angle  as  R  Is  varied.  Both  sets  of  Investigators 
adopted  a  common  fallacy  that  Interprets  dihedral  angles  In  terms  of  ratios 
of  y's  and  Ignores  the  torque  terms  that  result  from  the  N  dependence  of  y. 

Since  thermodynamics  requires  y  to  be  continuous  at  phase  transitions  no  discon¬ 
tinuity  In  the  ratio  of  y's  Is  possible.  As  a  result,  Gleiter's  Interpretation 
has  been  repeatedly  criticized  [2,10,11]. 

Much  of  the  discussion  has  centered  on  the  role  of  hypothetical  impurities.  But 
Impurities  are  chemical  components,  and  If  the  system  is  equilibrated  should  not 
affect  the  continuity  of  y.  Gleiter's  claim  that  the  angular  discontinuity  was 
reversible  Is  evidence  for  equilibrium. 

A  close  comparison  of  Gleiter's  Figures  4a  and  4b  [3]  reveal  evidence  for  the 
onet  of  facets  at  the  transition.  Coexistence  of  facets  Is  not  only  evidence  of 
phase  coexistence,  but  it  Is  also  evidence  of  torque  terms  resulting  from  the  N 
dependence  of  y.  As  can  be  readily  seen  from  £  vector  constructions,  there  is 
no  violation  of  thermodynamic  principles  by  a  report  of  a  discontinuous  change  In 
two  of  the  dihedral  angles  when  one  of  the  boundaries  at  a  trijunction  undergoes  a 
facetting  transition. 

General  Dissociation  Transitions  -  The  replacement  of  one  grain  boundary  with 
several,  whose  combined  misorientatlon  has  the  imposed  R,  and  whose  combined  y's 
Is  lower.  Is  a  special  case  of  a  more  general  dissociation  phenomenon.  When  the 
grain  boundary  dissociates  into  two,  a  grain  of  the  same  phase  has  formed  a  layer 
that  "wets"  the  original  grain  boundary,  lowering  Its  energy.  When  the  bismuth 
grain  melts  it  Is  replaced  by  two  Interfaces  plus  a  layer  of  a  different  phase 
which  wets  the  grain  boundary.  When  a  single  layer  wets  the  boundary  It  Is  easy  to 
understand  the  phenomenon  in  terms  of  thin  layers  of  three-dimensional  phases 
bounded  by  two  Interfaces  or  grain  boundaries  [41,42].  Often  the  layer  Is  not  a 
homogeneous  phase,  but  one  that  changes  continuously.  It  then  becomes  easier  to 
consider  the  boundary  a  region  of  continuous  transition  [43],  In  this  section  we 
consider  only  the  dissociation  through  the  means  of  a  wetting  layer  composed  of  a 
uniform  phase  and  two  Interfaces  or  boundaries. 

If  the  wetting  layer  Is  a  different  three-dimensional  phase,  we  expect  a  wetting 
transition  only  If  the  thermodynamic  cost  of  Inserting  this  phase  Is  negative. 
The'phase  Is  therefore  likely  to  be  one  that  is  not  too  far  from  being  stable  In 
bulk.  Since  grain  boundaries  can  not  exist  without  the  grains,  the  ultimate  grain 
boundary  phase  transition  occurs  when  the  grains  themselves  change  phase.  At  such 
a  phase  transition  some  grain  boundaries  will  be  wet  and  the  thickness  of  the 
wetting  layer  Is  without  limit.  If  conditions  are  changed  so  that  the  wetting 
phase  Is  slightly  unstable  In  bulk.  It  will  remain  at  these  boundaries  but  It's 
thickness  will  be  goverened  by  forces  between  the  two  Interface  boundaries. 
Eventually  as  the  wetting  phase  becomes  more  unstable,  it  may  disappear  from  the 
grain  boundary  In  a  phase  transition.  Grain  boundary  melting  transitions  have  been 
calculated  to  occur  hundreds  of  degrees  before  melting. 

In  addition  there  Is  a  singularity  at  the  melting  point  Itself.  Because  a  wetting 
layer  can  Increase  without  limit  at  the  bulk  phase  transition,  use  of  the  Gibbs 
adsorption  equation  Indicates  that  the  temperature  and  chemical  potential  coeffi¬ 
cients  of  y  must  diverge.  This  too  has  been  confirmed  In  a  grain  boundary  melting 
calculation  [40].  Wetting  transitions  are  not  confined  to  melting  and  dissoci¬ 
ation.  Solid  state  phase  changes  can  also  lead  to  wetting  layers  of  a  different 
solid  phase.  Surprising  results  [44]  on  grain  boundary  diffusion  In  Iron  might  be 
an  Indication  that  this  occurs. 

Wetting  occurs  universally  near  critical  consolutlon  points  because  y  of  the 
Interface  between  the  two  critical  phases  Is  extraordinarily  low  [42].  As  a  result 
since  the  grain  boundaries  of  one  of  the  critical  phases  will  have  a  lower  y, 
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grain  boundaries  of  the  other  critical  phase  will  be  wet.  As  two-phase  coexistence 
Is  approached  there  could  be  two  singularities;  first  either  an  abrupt  transition 
to  a  wetting  layer  or  a  discontinuous  change  in  adsorption,  followed  by  the  singu¬ 
larity  when  the  other  critical  phase  becomes  stable  in  bulk. 

Pi scussion  -  We  have  made  an  exact  analogy  between  the  orientation  of  a  grain 
boundary  normal  and  composition  in  a  three-component  system.  They  are  similarly 
constrained  by  experiment.  Only  their  averages  are  fixed  and  the  system  is  free  to 
lower  its  free  energy  by  facetting  or  forming  phases  with  different  compositions. 
The  same  convexity  rules  and  tangent  construction  hold  for  equilibria.  Phase 
diagrams  are  similar,  differing  mainly  in  phase  rule  modifications  because  of 
symmetry. 

These  considerations  led  to  an  obvious  definition  of  what  we  should  call  a  grain 
boundary  phase.  A  smoothly  curved  grain  boundary  is  like  an  inhomogeneous 
phase.  Phase  coexistence  occurs  usually  with  discontinuities  in  orientation. 

With  6+C  state  variables  for  grain  boundary  properties,  symmetry  has  an  Important 
role  in  bringing  the  observations  into  a  manageable  form.  It  also  plays  a  role  In 
helping  locate  interesting  features  in  phase  diagram,  in  sectioning  the  hyperspace 
and  In  predicting  deviations  from  the  phase  rule. 

Consistent  with  the  ternary  analogy,  phase  transitions  generally  occur  by  the 
appearance  of  new  orientations.  The  onset  of  facetting  is  a  common  observation. 
Often  it  Is  brought  on  by  the  addition  of  certain  chemical  components  [45].  This 
must  be  interpreted  as  a  phase  change. 

Dissociation  and  wetting  are  also  commonly  observed.  The  onset  of  wetting  Is  seen 
as  a  phase  transition  as  the  state  variables  are  varied.  It  has  been  found  as  R  Is 
varied  [25,38]  increasing  y  of  the  boundary  or  as  T  is  varied  [40,41]  changing  the 
stability  of  the  wetting  phase.  The  transition  from  dissociated  to  undissociated 
boundaries  has  not  been  pinpointed. 

Most  of  the  previous  discussions  of  grain  boundary  phase  transitions  have  impli¬ 
citly  assumed  congruent  phase  changes.  We  have  found  that  this  Is  a  rare  occur¬ 
rence  at  equilibrium.  Far  from  equilibrium  It  can  occur,  just  as  partitionless 
phase  changes  occur  In  the  alloy  analogue.  Most  of  the  expected  observations  in 
these  discussions  required  a  pinpointing  of  singularities  in  equilibrium  properties 
and  this  requires  observations  near  the  equilibrium  transition.  None  of  these 
methods  of  observation  are  both  easy  and  convincing.  By  comparison-  facetting  Is 
easy  to  observe,  and  so  closely  associated  with  grain  boundary  phase  transitions, 
that  we  have  chosen  to  define  phase  coexistence  In  terms  of  facets.  The  analogue 
to  nucleation  and  growth  and  to  spinodal  decomposition  has  been  made  [29]  in  a 
paper  that  recognized  the  equivalence  of  orientation  In  problems  of  facetting 
mechanism  with  those  of  composition  In  phse  changes. 

Phase  transitions  far  from  equilibrium  may  proceed  with  a  change  of  morphology. 

If,  however,  a  transition  Is  reported  to  occur  congruently  it  Is  Imperative 
that  the  transition  be  re-examined  close  to  equilibrium.  The  observation  may 
simply  have  been  a  relaxation  back  to  equilibrium  of  a  boundary  without  a  phase 
change.  At  high  temperature  boundaries  often  lose  the  adsorption  without  a  phase 
change.  Low  and  high  temperature  properties  often  seem  to  belong  to  different 
phases,  but  careful  examination  often  reveals  continuity  and  thermodynamic 
consistency  without  a  phase  change. 

In  a  polycrystalline  specimen  each  boundary  may  have  one  or  more  phase  transi¬ 
tions  but  these  will  occur  at  different  temperatures  for  different  boundaries. 
Observations  on  global  properties  of  polycrystalline  samples  will  smear  out  the 
phase  transitions  that  do  occur.  Micrographic  Investigation  for  facetting  changes 
seems  to  be  the  obvious  way  to  survey  all  the  many  boundaries  In  a  polycrystalline 
sample.  83 


The  most  important  conclusion  of  this  paper  is  the  clear  identification  between 
grain  boundary  phase  change  and  facetting.  It  has  been  incorporated  into  our 
definition  of  phase  change.  Reports  of  phase  changes  on  curved  boundaries  without 
the  facetting  evidence  are  not  credible.  The  experimental  test  required  by  this 
definition  is  simple  enough  that  we  should  now  know  when  we  have  a  phase  transition. 
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Fig.  4  :  Portions  of  an  N  diagram  straddling  a  mirror  plane.  Multi¬ 
phase  coexistence  is  marked  by  a  tie  triangle,  a  tie  trapezoid  and 
tie  lines.  Critical  points  between  mirror  related  phases  must  lie  on 
the  mirror. 
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